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ABSTRACT 


On the basis of wind speeds and accelerations derived from U.S. Navy sponsored 300-mb. constant level 
balloon or transosonde flights made during 1953, 1955, and 1956, statistics are presented on the magnitude of the 
ageostrophic wind and its variation with latitude and wind speed. These statistics indicate that at 300 mb. 
the average angle between wind and geostrophic wind is 11 degrees and the mean magnitude of the vector deviation 
between wind and geostrophic wind is 12 m. sec.“. The data also show that, through the use of the geostrophic 
and gradient wind approximations, half the time errors greater than 29 percent and 11 percent, respectively, are 


introduced into the derived results. 


1. INTRODUCTION 


Previous evaluations of the magnitude of the ageo- 
trophic wind have in large part been based on comparisons 
of wind data with geostrophic wind data obtained from 
yhoptic weather maps, as illustrated by the works of 
Neiburger et al. [1] and Machta [2]. A few evaluations 
of the ageostrophie wind have been based on estimations 
ifthe partial wind derivatives in the equations of motion, 
willustrated by the work of Godson [3]. Durst and Gil- 
et [4] and Neiburger and Angell [5] evaluated the 
igeostrophic wind by determining the value of the indi- 
dual wind derivatives from constant level balloon flights 
ul introducing this value into the equation of motion. 
‘ince limited constant level balloon data were available 
ranalysis in the latter two papers, it seemed desirable 
“carry through further analysis upon the receipt of more 


“stant level balloon or transosonde data, This paper 
en 
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presents statistics on the magnitude of the ageostrophic 
wind, the variation of this magnitude with latitude and 
speed, and indicates the errors made in assuming the wind 
to be geostrophic or gradient, all based on 635 evaluations 
of the ageostrophic wind from 28 transosonde flights made 
at 300 mb. Of these 28 transosonde flights, 8 were 
launched from Minneapolis, Minn. in 1953, 4 were 
launched from Vernalis, Calif. in 1955, and 16 were 
launched from Oppama, Japan in 1956. The average 
duration of these flights was 65 hours and the average 
trajectory length was 4,700 nautical miles. The flights 
were positioned at 1-hourly or 2-hourly intervals by 
means of the excellent Federal Communications Commis- 
sion (FCC) radio-direction-finding network. Anderson 
[6] found, by comparison of FCC positions with transo- 
sonde positions determined by photographs of the under- 
lying terrain, that within the United States two-thirds of 
the time the FCC positions are in error by less than 20 
nautical miles. Over the Atlantic and Pacific Oceans one 
would expect this error to be larger owing to the greater 
distances between tracking stations and balloons. 
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2. DETERMINATION OF THE WIND 
AND AGEOSTROPHIC WIND 


With minor approximations, the vector equation of 
motion may be written : 


on =f(V—V,) Xke-+F (1) 


where V is the horizontal component of the velocity, V, 
is the geostrophic velocity, f is the Coriolis parameter, F 
is the frictional force per unit mass, and k is the unit 
vertical vector. The term on the left hand side of (1) 
represents the change of horizontal velocity following 
the three-dimensional path of an air parcel and may be 
expanded into three terms involving partial derivatives, 


dV_aVv av 


where w is the vertical velocity and V is the horizontal 
differential operator. The constant level balloon data 
determine the sum-value of the first two terms on the 
right hand side of (2). The sum of these two terms (lo- 
cal change of velocity plus horizontal advection of veloc- 
ity) is hereafter called the partial acceleration. The 
value of the third term on the right hand side of (2). 
the vertical advection of velocity, may be approximated 
along a constant level balloon trajectory by determining 
the vertical motion by the adiabatic technique (vertical 
motion proportional to the change of temperature follow- 
ing the balloon and inversely proportional to the devia- 
tion of the lapse rate from the process lapse rate along 
the trajectory), and by determining the vertical wind 
shear from rawin data. Based on a limited number of 
computations from the 1953 flights and reported in refer- 
ence |5], the magnitude of the vertical advection of veloc- 
ity averages about 25 percent of the magnitude of the 
partial acceleration, and thus is certainly not negligible. 
Nevertheless, in this paper it is assumed that the partial 
acceleration data derived from the transosondes repre- 
sents a reasonable approximation to the total accelera- 
tion since the bulk of the transosonde positions are over 
oceans where it is impractical to estimate the vertical 
advection of velocity by the techniques described above. 

In order to pass from the value of the acceleration to 
the value of the geostrophic deviation it is necessary, 
as seen from (1), to neglect friction. While attempts 
have been made to estimate the magnitude of the fric- 
tional force from a combination of radiosonde and transo- 
sonde data, the evaluations, particularly over the oceans, 
are far too crude for inclusion here. Hence, an additional 
degree of uncertainty (in addition to the neglect of the 
vertical advection of velocity) enters into the values of 
the geostrophic deviations presented below. 

With the neglect of friction and the vertical advection 
of velocity we find, after taking the dot product of (1) 
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with the unit vector along and normal to the y ind, that, 
respectively, 


dV/dt=fV, sin i (3) 
Vd6/dt= —f(V—V, cos 1) (4) 


and 


where d6/dt is the angular velocity of an air parcel and; 
is the angle of indraft (the angle between wind ani 
geostrophic wind). Taking into account the findings of 
Neiburger and Angell with respect to the time interya 
for which the velocity and acceleration should be coy. 
puted in order to avoid large errors in these parameters 
due to radio-direction-finding errors in the positioning of 
the transosondes, the speed was determined from a 4-how 
transosonde displacement; the tangential acceleration 
(dV /dt) was determined from the change in 4-hour aver. 
age transosonde speed in 8 hours; and the normal x. 
celeration (Vd6/dt) was determined as the product of the 
mean transosonde speed during an 8-hour time interval 
and the change of 4-hour average transosonde direction in 
8 hours. 

Upon eliminating the geostrophic wind speed between 
(3) and (4), the angle of indraft can be expressed in 
terms of three known variables; namely, the tangential a. 
celeration, the normal acceleration, and the Coriolis x. 
celeration (fV). From knowledge of the angle of in. 
draft, the geostrophie wind can be evaluated from eithe 
(3) or (4). The cross-contour component of the ageo 
strophic wind (V sin 7) can then be determined as wel 
as the along-contour component of the ageostrophie win 
(V cos i—V,). In the diagrams in this paper V—J, is 
substituted for V cos i—V, since, for the magnitudes of 
the angles of indraft obtained at 300 mb., the two quant: 
ties are practically identical. 


3. VELOCITY AND?AGEOSTROPHIC 
VELOCITY STATISTICS 


In this section and in sections 4 and 5 we show bj 
means of histograms and cumulative frequency curve 
the magnitudes of various velocity and ageostorophic ve 
locity parameters obtained from the transosonde flights 
Where comparison of two parameters is desirable, histo 
grams and cumulative frequency curves are superimpos’ 
on the same diagram. 

The magnitude of the 4-hour average wind speed # 
300 mb. for this series of transosonde flights is indicated 
in the upper diagram of figure 1. Based upon 755 calet: 
lations, the mean wind speed is 38 m. sec.~? while the mode 
is 20-30 m. sec.. Ten percent of the winds exceed 67 ™ 
sec. The maximum wind speed determined from 4+ 
hour transosonde displacement is 150 m. sec’. This value 
occurred on Flight 23 approximately 1,000 miles eas 
southeast of Tokyo in January 1956, and is undoubted! 
one of the largest, reasonably reliable, values of wind 
speed ever recorded. 
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-Absolute magnitude (top) and algebraic magnitude 
(bottom) of angle of indraft, 7, at 300 mb. 


FIGURE 2.- 





The magnitudes of the zonal and meridional wind com- 
ponents are indicated in the lower diagram of figure 1. 
These data demonstrate the well known fact that the mean 
value of the zonal wind component (here 30 m. sec.-*) is 
considerably larger than the mean value of the meridional 
wind component (17 m. sec.*). 

The absolute magnitude of the angle of indraft, indi- 
vated in the top diagram of figure 2, is based on 635 caleu- 
lations. The mode is indistinct and lies in the range 0-6 
degrees, while the median is 8 degrees and the mean is 11 
legrees. In 15 percent of the cases the angle exceeds 20 
(legrees. ‘The considerable difference between median and 
lean is attributable to a relatively few large values of the 
ingle of indraft of doubtful validity. For comparison, 
todson found by substituting synoptic wind data at 700 
mb. in the equation of motion a mean value of 14 degrees 
for the angle of indraft. The larger mean mean value 
for his data would be expected since the magnitude of 
the angle of indraft is inversely proportional to the wind 
‘peed, as shown in section 6 of this paper. 

In the lower diagram of figure 2 is shown the algebraic 
magnitude of the angle of indraft, considered positive 
when the ‘low is toward low pressure and negative when 


the flow is toward high pressure. It is evident from the 
cumulative frequency curves in this diagram that in the 
mean the angle of indraft is of greater magnitude when 
the flow is toward high pressure than when the flow is 
toward low pressure. It is believed that this tendency is 
at least partially connected with the occasional loss of 
geostrophic control in regions of inertial instability in the 
deceleration region downstream from sharply curved 
ridges. 

The magnitude of the vector geostrophic deviation 

V-—V, | is indicated in the upper diagram of figure 3. 
The mode is 2-4 m. sec.’, the median is 10 m. see.' and 
In 15 percent of the cases the 
For comparison, 


the mean is 12 m. sec.’ 
deviation is greater than 20 m. sec". 
Machta found an average value of 13 m. sec.’ for the vector 
geostrophic deviation at 300 mb, during the winter 
months. His data were obtained by comparing the wind 
and pressure gradient (geostrophic wind) on analyzed 
synoptic maps. 

The absolute magnitudes of the cross-contour (J" sin @) 
and along-contour (V—TJ’,) components of the vector 
geostrophic deviation are given in the lower diagram of 
figure 3. In the case of the cross-contour component, the 
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Fiacure 3.— Magnitude of the 300-mb. vector geostrophic deviation, 
|\V—V,|, (top) and its cross-contour, |V sin i], and along-contour, 
|V—V,|, components (bottom). 
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Fiaure 4.—Algebraic magnitude of the 300-mb. cross-contour flow, 
V sin 7, (top) and along-contour component of the vector geo- 
strophic deviation, V—V,, (bottom). 





mode, the median, and the mean are respectively 0-4 m. 
sec.", 4 m. sec.-?, and 6 m. sec.“ whereas for the along- 
contour component they are, respectively, 0-4, m. sec.”’, 
7 m. sec.-', and 9 m. sec.-*. In the case of the cross-con- 
tour component, 10 percent of the values exceed 13 m. sec.~’ 
whereas in the case of the along-contour component, 10 
percent of the values exceed 19 m. sec.*. For comparison, 
Neiburger et al. found at 700 mb. a mean value of 3 m. 
sec." for the along-contour component of the vector 
geostrophic deviation. This smaller value would be ex- 
pected since the magnitude of this component is a func- 
tion of the wind speed, as will be shown below. 

The algebraic magnitude of the cross-contour flow is 
shown in the upper diagram of figure 4. In agreement 
with the findings for the angle of indraft, the cumulative 
frequency curves indicate that the magnitude of the cross- 
contour flow is greater when the flow is toward high pres- 
sure than when the flow is toward low pressure. 

The algebraic magnitude of the along-contour com- 
ponent of the vector geostrophic deviation is given in the 
lower diagram of figure 4. In this diagram the differ- 
ence between the cumulative frequency curves is striking 


indeed, with the difference between the wind speed and 
geostrophic wind speed being much larger, in the mean, 
when the flow is cyclonic than when the flow is anti- 
cyclonic. The respective median values are —8 m. sec.’ 
and 6 m. sec.-". It may be that this asymmetry is due to 
the existence of a type of upper bound on the degree to 
which the wind may be supergeostrophic on ridges (twice 
the geostrophic value), whereas there is no such bound in 
troughs. 


4. ERRORS RESULTING FROM THE USE OF 
GEOSTROPHIC AND GRADIENT WIND EQUATIONS 


The magnitude of the partial acceleration 
(2V/2t+V-VV) in comparison with the Coriolis ae 
celeration (fV) is of interest since the geostrophic %& 
sumption states that the partial acceleration is negligible 
in comparison with the Coriolis acceleration. The mag 
nitude of the ratio of partial acceleration and Coriolis 
acceleration is given in the upper diagram of figure > 
The mode of the ratio is roughly 0.10-0.20 while the 
median is 0.29 and in 10 percent of the cases the rato 
exceeds 0.83. In other words these data suggest that, 
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Figure 7.—Algebraic ratio of the angular velocity of the wind and 
the angular velocity of the earth about the local vertical (d@/dt)/Q, 
at 300 mb. 





through the introduction of the geostrophic approxima- 
tion, half the time an error exceeding 29 percent is intro- 
duced into the derived results. 

The magnitudes of the ratios of zonal and meridional 
components of partial acceleration and Coriolis accelera- 
tion are given in the lower diagram of figure 5. These 
data suggest that the neglect of the acceleration in the 
wnal component of the equation of motion results in an 
wtenable approximation, since for this component the 
node of the ratio is 0.30-0.40 and the median is 0.41. In 
the case of the meridional component of the equation of 
motion, however, the ratio of partial acceleration and 
Voriolis acceleration has a mode of only 0-0.10 and a 
iedian value of only 0.25. In this diagram the relatively 
large percentage of cases where the ratio exceeds 1.00 is 
ue mainly to the existence of “cusp” points in the tra- 
ectories. At such points the velocity component may go 
toward zero at a time when the component of the partial 
weeleration is far from zero. 

The ms enitude of the ratio of tangential and normal 
‘omponer s of partial acceleration is of interest since the 


gradient wind approximation assumes the tangential com- 
ponent to be negligible in comparison with the normal 
component. It is seen from figure 6 that the mode of this 
ratio is 0-0.20 while the median is 0.67. In 38 percent 
of the cases the tangential component of the partial accel- 
eration is larger than the normal component. From con- 
sideration of the median ratio of the tangential and nor- 
mal components of partial acceleration and the median 
ratio of partial acceleration and Coriolis acceleration we 
estimate that through the use of the gradient wind ap- 
proximation half the time an error exceeding 11 percent is 
introduced into the derived results. 


5. FREQUENCY OF OCCURRENCE OF ‘‘ABNORMAL” 
FLOW 


In figure 7 is indicated the ratio of the algebraic angular 
velocity of the air parcel or transosonde and the angular 
velocity of the earth about the local vertical. It is seen that 
in 27 percent of the cases the angular velocity of the air 
parcel exceeds that of the earth about the local vertical 
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left), and of the vector geostrophic deviation, |W—V,|, (upper 
left) and its along-contour, |V—V,|, (upper right) and cross- 
contour, |V sin 7|, (lower right) components at 300 mb. 


and of the vector geostrophic deviation, |\W—V,,|, (upper left) and 
its along-contour, |V—V,|, (upper right) and _ cross-contour, 
|V sin t|, (lower right) components at 300 mb. 





when the angular velocity is positive (cyclonic flow), but 
that this criterion is satisfied in only 9 percent of the 398 
cases when the angular velocity of the air parcel is nega- 
tive (anticyclonic flow). Thus in 5 percent of the total 
635 cases the flow allegedly possesses anticyclonic rotation 
in space. Flow which satisfies this criterion is called “ab- 
normal” or “anomalous” flow by Holmboe et al. [7] and 
Gustafson [8]. 


6. VARIATION OF AGEOSTROPHIC PARAMETERS 
WITH LATITUDE AND SPEED 


In this section and in section 7 the data are presented in 
the form of group means, regression lines, and correlation 
coefficients (figs. 8-11). On either side of the group means, 
lines have been drawn extending a distance equal to two 
standard errors from the mean. If the data are randomly 
drawn from a normal population there is only a 5 percent 
chance that the true group mean lies outside the extent of 
these lines. An estimate of significance of variability at 
the 5 percent level can then be obtained by noting whether 
a straight vertical line can be drawn so as to intersect all 
these lines extending either side of the group mean. The 
regression lines and correlation coefficients serve to yield 


an indication of the overall trend but it is not meant to 
imply thereby that these trends are actually linear in 
nature. It is also possible to determine the correlation 
coefficient which would indicate significance at the 5 per- 
cent level. Based on the 2 transformation presented in 
Hoel [9], it is found that for the number of cases available 
here a correlation greater than +0.08 is significant at the 
5 percent level. 

With regard to the variation of the parameters with 
respect to latitude it must be emphasized that the results 
are biased since the transosondes sample chiefly cyclonic 
flow patterns at southerly latitudes and anticyclonic flow 
patterns at northerly latitudes. If transosondes were Ie 
leased from stations along a given meridian but at differ- 
ent latitudes, it is probable that different mean values 
would be found from those presented below. 

Figure 8 gives the variation with latitude of the ab- 
solute magnitude of the angle of indraft and the vector 
geostrophic deviation and its along-contour and cross-col 
tour components. All four parameters decrease in mag 
nitude with increasing latitude, with the smallest decreas 
being found for V sin ¢ (r=—0.12). The group meals 


show that the variability is ragged. The magnitude of 
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the vecto. geostrophic deviation varies from a mean of 
{) m. sec’ at northerly latitudes to 14 m. sec. at south- 
aly latitudes with most of the change occurring abruptly 
at latitude 35°. The absolute magnitude of the angle of 
indraft varies from a mean value of about 8 degrees at 
jortherly latitudes to a mean value of 13 degrees at south- 
erly latitudes, while the absolute magnitude of V sin 7 
varies from a mean value of 45 m. sec.~? in northerly 
latitudes to a mean value of 7-8 m. sec. in southerly 
latitudes. 

Figure 9 gives the variation with speed of the absolute 
magnitude of the angle of indraft and the vector geo- 
trophic deviation and its along-contour and cross-contour 
components. While the angle of indraft shows a signifi- 
cant decrease in magnitude with increase in speed 
(r=-0.27) the ageostrophic parameters show significant 
increases in magnitude with increase in speed. It is of 
interest to note that the magnitude of the cross-contour 
fow increases with increase in speed despite the counter- 
ating tendency of the angle of indraft. The vector 
gostrophic deviation varies in average magnitude from 
$m. sec.? at low speeds to 18 m. sec.‘ at high speeds. The 
absolute magnitude of the angle of indraft varies from 17 
legrees at low speeds to 7 degrees at high speeds, while 
the absolute magnitude of the cross-contour flow varies 
from 4m. sec.~? at low speeds to 11 m. sec.-' at high speeds. 


7. VARIATION WITH LATITUDE AND SPEED 
OF ERRORS RESULTING FROM THE USE OF 
GEOSTROPHIC AND GRADIENT WIND EQUATIONS 


The two diagrams in figure 10 give the variations with 
latitude of the ratio of partial acceleration and Coriolis 
wceleration and the ratio of tangential and normal com- 
ponents of partial acceleration. As would be expected, the 
ratio of partial acceleration and Coriolis acceleration 
varies greatly with latitude (r—=— 0.32), ranging from a 
ean Value of 0.30 at northerly latitudes to a mean value 
of 0.70 at southerly latitudes. The group means indicate 
that the ratio is almost constant north of latitude 40°, 
with a rapid, nearly linear increase in the value of the 
ratio as one progresses southward from latitude 40°. The 
ritio of the tangential and normal components of the 
partial acceleration shows little variation with latitude 
(r=—0.03). The group means also indicate no significant 
variations except at latitude 60° where the large values of 
thenorma! components of partial acceleration on the crests 
of trajectories overbalance the tangential components. 

The two diagrams in figure 11 give the variations with 
speed of the ratio of partial acceleration and Coriolis ac- 
‘sleration and the ratio of tangential and normal com- 
ponents of partial acceleration. As would be expected, 
the ratio of partial acceleration and Coriolis acceleration 
shows a significant decrease in magnitude with increase 
N Speed (7=—0.38) with the mean value of the ratio 
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varying from 0.66 at low speeds to 0.20 at high speeds. 
The ratio of tangential and normal components of partial 
acceleration shows no significant variation with speed 
(r= —0.02), the group means at all speeds resting close to 
0.70. 


8. CONCLUSION 


Statistics based on 635 evaluations from 300-mb. con- 
stant level balloon data indicate that through the use of 
the geostrophic wind approximation half the time an er- 
ror greater than 29 percent will be introduced into the 
derived results, whereas through the use of the gradient 
wind approximation half the time an error greater than 11 
percent will be introduced into the derived results. On 
a percentage basis both the geostrophic and gradient wind 
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approximations are about twice as bad at latitude 20° as 
at latitude 60° and about twice as bad at speeds of 10 m. 
sec. as at speeds of 50 m. sec... The magnitudes of 
angles of indraft and ageostrophic parameters are in fair 
agreement with results obtained by other techniques, with 
the transosonde data yielding a mean angle of indraft of 
11 degrees, a mean cross-contour flow of 6 m. sec.-', a mean 
deviation between wind and geostrophic wind speed of 
9 m. sec., and a mean vector geostrophic deviation of 
12 m. sec."’. 


REFERENCES 


1. M. Neiburger et al., “On the Computation of Wind from Pres- 
sure Data,” Journal of Meteorology, vol. 5, No. 3, June 1948, 
pp. 87-02. 

2. L. Machta, A Study of the Observed Deviations from the 
Geostrophic Wind, ScD thesis, Dept. of Meteorology, Mas- 
sachusetts Institute of Technology, 1948. (Unpublished). 


3. W. L. Godson, “A Study of the Deviations of Wind Speeds and 


Directions from Geostrophic Values,” Quarterly Jowrng) ,, 
the Royal Meteorological Society, vol. 76, No. 327, Jan, 1% 
pp. 1-15. 


. C. L. Durst and G. H. Gilbert, “Constant-Height Balloons 


Calculation of Geostrophic Departures,” Quarterly Joury,, 
of the Royal Meteorological Society, vol. 76, No. 327, Jy 
1950, pp. 75-86. 


. M. Neiburger and J. K. Angell, “Meteorological Applications y 


Constant-Pressure Balloon Trajectories,” Journal of Metey,. 
ology, Vol. 13, No. 2, Apr. 1956, pp. 166-194. 


3. A. D. Anderson, “A Study of the Accuracy of Winds Deriyy 


from Transosonde Data,” (Memo. Rep. 498), Naval Resear) 
Laboratory, Washington, D.C., 1955, 10 pp. 


. J. Holmboe, G. E. Forsythe, and W. Gustin, Dynamic Meteor. 


ology, New York, John Wiley and Sons, 1945, 378 pp. 


._ A. F. Gustafson, “On Anomalous Winds in the Free Atmm 


phere,” Bulletin of the American Meteorological Society, yo 
34, No. 5, May 1953, pp. 196-201. 


. P. G. Hoel, Introduction to Mathematical Statistics, New Yor 


John Wiley and Sons, 1947, 258 pp. 








AY 19) 


pds ani 
rnal of 
LD. 195) 


loons— 
Tourna| 
7, Jar, 


ions of 
VW eteoy - 


derive 
search 


ete Or. 


Atmos. 
ty, Vol. 


York, 





May 1969 





MONTHLY WEATHER REVIEW 171 


ON BAROTROPIC AND BAROCLINIC MODELS, WITH SPECIAL EMPHASIS 
ON ULTRA-LONG WAVES 


A. WIIN-NIELSEN, Air Weather Service 
Joint Numerical Weather Prediction Unit, Suitland, Md. 
(Manuscript received April 1, 1959; revised May 22, 1959] 


ABSTRACT 


The problem of control of the ultra-long waves in numerical prediction is considered. 
It is shown that the vertical variation of the horizontal wind and the static 


discussion of the one-level forecasts. 


stability are the main factors in determining the value of divergence at 500 mb. 


Sections 2-4 contain a 


An independent estimate of the 


size of the term governing the ultra-long waves in the atmosphere agrees well with the one determined by 


Cressman on an empirical basis. 


Section 5 points out that any two-parameter model has to contain an effect similar to the one contained in the 
one-parameter model controlling the ultra-long waves. A modification of a two-parameter model is made in such 


a way that the ultra-long waves are controlled. 


Section 6 describes a perturbation analysis of the model developed in Section 5 in order to investigate the effect 


of the modification also on the shorter waves. 


It is found that a certain stabilization of the shorter waves 


is produced. Baroclinic instability and growth rate are investigated. 
Sections 7 and 8 contain a justification of certain approximations used in the earlier sections regarding the 


vertical variation of static stability and the profile of vertical velocity. 


1. INTRODUCTION 


The Joint Numerical Weather Prediction Unit has re- 
cently improved the 500-mb. barotropic forecasts by 
paying special attention to the behavior of the very long 
atmospheric waves. The non-divergent barotropic fore- 
casts were first empirically corrected by Wolff [18], who 
ued the observed fact that the very long atmospheric 
waves do not change either position or ampltitude to any 
great extent from day to day. He therefore forced the 
waves with wave-number one, two, and three to be sta- 
tionary throughout the forecast period. Shortly after, 
Cressman [6] used a special case of Phillips’ [13] two- 
layer model to show that a model atmosphere consisting 
of two homogeneous layers, where the motion in the lower 
fluid is parallel to the contour of the interface, and where 
the motion in the upper fluid is negligible, would decrease 
the retrogression of the very long atmospheric waves to a 
considerable extent. The forecast equation used by Cress- 
man is the same as the one used by Bolin [2] in his so- 
called tropopause model. Both of these authors had dif- 
ficulties in the determination of the proper value for 
vertain constants appearing in the prognostic equation, 
itinly because a model atmosphere consisting of two 
homogeneous layers is rather difficult to “translate” to the 
real atmosphere. Cressman decided therefore to deter- 
mine the value of the coefficient by computing a number 
of forees sts on the same data with different values of the 
‘oeflicien: and then choosing the coefficient that gave the 
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best forecast. This procedure may not lead to the proper 
value of the coefficient, because a minimization of forecast 
height errors with respect to the coefficient might tend to 
compensate for correlated, but physically unrelated er- 
rors, as was also mentioned by Cressman. 

In view of the situation mentioned above it seems to be 
in order to try to formulate the barotropic, divergent 
model in such a way that the value of the constants can 
more easily be related to the real atmosphere. Sections 
2 through 4 of this paper are devoted to this purpose. 

An analysis of baroclinic instability and especially of 
the phase-speed of waves in a baroclinic atmosphere, as 
for instance made by Eliassen [9], shows that the phase- 
speed for large values of the wavelength approaches the 
Rossby speed for a non-divergent atmosphere. For large 
values of the wavelength one actually obtains two solu- 
tions for the phase-speed. One of them is not likely to 
be observed in the atmosphere because it corresponds to a 
situation where the temperature wave is out of phase with 
the pressure, corresponding to warm troughs and cold 
ridges. The other solution is the one approaching the 
Rossby speed. 

This result shows that one will find about the same dif- 
ficulties as in the non-divergent barotropic model with 
respect to the ultra-long waves, when one integrates a 
baroclinic model of the usual type over an almost hemi- 
spheric region. It is therefore a necessity to correct the 
baroclinic forecast equations in such a way that the 
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retrogression of the ultra-long waves is reduced to a con- 
siderable extent. Section 5 describes a simple two- 
parameter model where this effect is included. 

Recently Burger [3] pointed out that the vorticity 
equation loses its prognostic value when the wavelength is 
of the same order of magnitude as the radius of the earth. 
As long as we base our forecasts on the vorticity equation 
we can therefore probably not do anything better than 
express the quasi-stationary character of the ultra-long 
waves. The two models presented here do therefore not 
claim to have any skill in the regime of the ultra-long 
waves. The problem which has been attacked is to in- 
corporate effects in barotropic and baroclinic models 
which control the long waves and at the same time to in- 
vestigate the changes which may be caused by such effects 
on the shorter waves. 


2. DERIVATION OF THE PROGNOSTIC EQUATION 
IN THE BAROTROPIC CASE 


In the derivation we shall consider the vorticity equa- 
tion in the form 


SEV vate $e (2.1) 
where V is the horizontal wind vector, { is the vertical 
component of the relative vorticity, »>=¢+/f, f is the Cori- 
olis parameter, and o=dp/dt the vertical velocity in 
a coordinate system with pressure as the vertical co- 
ordinate. 

This form of the vorticity equation, where the value of 
the Coriolis parameter in the divergence term is a stand- 
ard value, f==/,, is consistent with certain general integral 
constraints, as shown by the author [17], when the hori- 
zontal wind in the vorticity advection term is non- 
divergent. The derivation of the barotropic vorticity 
equation is usually based upon the equivalent-barotropic 
atmosphere (Charney [5]), where the horizontal wind 
in the complete atmosphere is assumed to vary vertically 
in strength, but not in direction. It is then shown that 
the simple barotropic vorticity equation can be applied 
to the vertical mean flow in the atmosphere provided the 
effect of surface pressure changes is neglected. The 
assumption 

V=A(p)Vv* (2.2) 
which defines the equivalent-barotropic atmosphere and 
where V* for instance may be taken as the 500-mb. flow, 
will in general apply with good accuracy in a layer 
around the 500-mb. level, while the accuracy becomes 
smaller the farther we are from the 500-mb. level. We 
shall in the following make use of (2.2) only in a thin 
layer around the 500-mb. level. Equation (2.2) expresses 
then, according to the geostrophic thermal wind relation, 
that the thermal wind is parallel to the wind itself, or in 
other words that the isotherms at 500 mb. are parallel to 
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the contour lines. The last requirement is probubly fy}. 
filled to the greatest extent in the mid-troposphere, whi 
temperature advection by the horizontal wind is greater 
both higher up and lower down in the troposphere, ge. 
cording to observations. 

The equivalent-barotropic atmosphere is usually treated 
without any reference to the thermodynamics of the 
atmosphere. If, however, the relation (2.2) is approxi. 
mately satisfied at least around the 500-mb. level! it muy 
mean that the local temperature changes to the larges 
extent are compensated by vertical velocities, and that the 
adiabatic equation therefore takes the form 


5 (gp) tee=0, o=—a dp (2.3) 


at the 500-mb. level. In (2.3) ¢=gz is the geopotential, 
a the specific volume, and 6 the potential temperature, 
The expression (2.2) applies actually to the distribution 
of the horizontal wind. If the horizontal wind is assumed 
quasi-geostrophic, we obtain from (2.2): 
0'/2¢\_dA do" 
ot (3 dp Ot 


olnée 





(2.4 


The vertical velocity can then be obtained from (23 
in the form: 

— —¢Aldp a9" 

= “ ry, (2.5 


We shall further in the thermodynamic equation (2.5) 
replace the geopotential ¢* by a stream function /’ 
satisfying the relation 


¢* 
+", (2.6) 
wry, 


The implication of this assumption in the thermody- 
namic equation, and only here, has been discussed by 
Phillips [15]. Introducing (2.6) in (2.5) we obtain: 

x* 
dA/dp fo oy* (27) 


o ot 





In order to obtain the prognostic equation from (2.1) 
it is now necessary to evaluate 0w/dp from (2.7). Letus 
assume for simplicity that dA/dp is a constant. The main 
problem is then in which way the vertical stability o varies 
with pressure. It will be shown (section 7) that a varit- 
tion of the vertical stability expressed by the formula 


-2 
o=ot (B) » P=500 mb. (2.8 


describes the vertical variation of temperature in a meal 
atmosphere with good accuracy. In (2.8) o* is the state 
stability at 500 mb. The accuracy is especially gout 
around 500 mb., because the vertical stability at 500 mb. 
is used to fit the distribution (2.8) to the atmosphere. 

From (2.7) and (2.8) we obtain then, provided 4(?) 
varies linearly with pressure : 
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=) ___ 2fo 8 (3)% oy* (2.9) 
which, inserted in (2.1), leads to the prognostic equation: 
S (ryt) + JV", n*)=0 (2.10) 
where 2f2 ae 
_ 2fo 
=P ($s P (2.11) 


3, DISCUSSION OF THE PROGNOSTIC EQUATION 


The Helmholz term appearing in (2.10) is in this for- 
mulation due to two factors, the increase of the horizontal 
wind with height and the decrease of temperature with 
height in the atmosphere. Both of these factors are well 
known from observations, and it is therefore somewhat 
easier to determine the proper value of g to be used here 
than in Bolin’s and Cressman’s prognostic equations. In 
order to determine g we need a value for o* and for 
(-dA/dp). 

The value for o* may be determined from the expression 


Olnée RT 


Sp =opF OM —y) (3.1) 


c= --a2a => 


where 7? is the gas constant, ya the dry adiabatic lapse rate, 
and y the local lapse rate. 

It has been the general experience among people work- 
ing With models in numerical prediction that the standard 
atmosphere is too unstable. Using y=0.5° C. per 100 m., 
which corresponds to about 80 percent of y in the standard 
atmosphere, we obtain : 


=4.2 MTS units. (3.2) 


We have already assumed that A(p) varies linearly 
with pressure in the derivation of (2.9). We may now 
obtain dA /dp by noting that A (p*)—1 according to (2.2) 
und further assume that A(2P) =0.2, which means that 
the surface wind is about 20 percent of the 500-mb. wind, 
ai assumption used, for instance, in the incorporation of 
the mountain effect in the barotropic model. This gives: 


—dA/dp=0.8/P=16 x 10-* (3. 3) 
and 
q@i.5X10-” (3. 4) 


The prognostic equation (2.10) may now be compared 


with Cressman’s prognostic equation reproduced here 
without the term accounting for the mountain effect : 


(vv rg v)+ J(¥*,n*) =0 (3.5) 


In this « uation the symbols have the same meaning as 
lefore. “J * is a representative mean value of the stream 


function obtained through the solution of the balance 
equation, while » is the coefficient, which actually is equal 


, 
tol fo-£) where / is the height of the interface, z the 


height of the 500-mb. surface, and p’ and p the densities of 
the upper and lower fluid. We are now in the position to 
evaluate the value of », which would correspond to the 
estimate of g given in (3.4), if we linearize the coefficient 
in (3.5). We obtain then: 


_gz* 
gee 


This independent estimate of » is about equal to the 
highest value of » used by Cressman in his empirical 
determination of the coefficient. Cressman varied » from 
0 to 8. It will be seen that » is inversely proportional to 
o*. A value of o* higher than the one used here cor- 
responding to y=0.5°C./100 m. is probably not possible 
as a mean value for the atmosphere. If we therefore 
accept this value of o*, it is found that values of , less 
than 8 can be obtained only with —dA/dp smaller than 
the value used here. It is interesting to note that a value 
of »=4, which is the value used at the moment in the 
operational forecasting, would correspond to —dA/dp 
=(0.4/P. This would mean that the surface wind should 
be about 60 percent of the 500-mb. wind, which sounds 
unreasonable. On the basis of this discussion one is 
therefore tempted to conclude that a value of » somewhat 
greater than 4 ought to be used. 


q=~8 (3.6) 


4. A SIMPLE WAVE ANALYSIS 


It is of interest to make a simple wave analysis of the 
prognostic equation (2.10) in order to obtain a first idea 
of the influence which different values of g will have on 
the motion of long and short waves in the atmosphere. 

Consider a stream function 


¥(2,y,t) = —Uy+ Ae*®?- (4.1) 


where the superscript * now has been dropped and where 
U’=constant is the zonal component of the wind. This 
stream function will satisfy (2.10) provided the phase 
speed c satisfies the relation 


— Bik? _ Cnp (4 9 ) 
= ae i+q/e . 





In the derivation of (4.2)8=df/dy has been considered as 
constant and is evaluated at 45° N. k=2x/L, where L is 
the wavelength. If N is the number of waves around the 
hemisphere, i.e., 2xR cos g=N- L, we may also write 
(4.2) in the form 


U—29R cos y/N? 


m= gh con gINT (4.3) 











a ee 
tt ee ST aT. oe cP 











Samet 
012 8 








+ a 5 10 4 2 a, 


Fieurgk 1.—VPhase velocity as a function of number of waves 
around the hemisphere. The solid line corresponds to g=0, the 
dashed line to gq=(3/4)10-" and the dotted line to gq= (3/2)10™. 


The values of ¢ shown in figure 1 in m. sec.-? are com- 
puted for /=20 m. sec.", p=45° N., and for g=0, 
(3/4) X 10°, and (3/2) X10-*, which would correspond 
to approximately »=0, 4, and 8. It is seen that 
qg= (3/2) X 107, the value estimated in this study, reduces 
the retrogression further than the operational value. For 
wave number 1 the phase-velocity changes from about 
—19 m. sec.-' to about —11 m. sec.-?, while on the other 
hand ¢ is changed from 12 m. sec.* to 10 m. sec.-* for 
N=8, corresponding to a wavelength of about 3500 km. 
It is therefore likely that a change of g from the present 
value, (3/4) X10", to the value determined here, 
(3/2) X10-*, would not change the forecast appreciably, 
except for the ultra-long waves, which in fact are 
better forecasted by the greater value of g, according to 
Cressman. 


5. THE BAROCLINIC CASE 


The incorporation of a control of the ultra-long waves 
in a baroclinic model is complicated by the fact that we 
want the corrected model to keep essentially the same 
properties it had for the shorter waves. A pronounced 
difference between a barotropic and a baroclinic model is 
that. the baroclinic model contains a mechanism for in- 
stability by which waves may grow in amplitude. Once 


we have decided in which way we want to incorporate the 


MONTHLY WEATHER REVIEW 


May 19% 


control of the long waves we have therefore to in vestigate 
whether the phase-speed and the instability criteria for the 
modified model are left essentially unchanged, or if they 
are changed, we have to investigate whether the « change j is 
in a direction which is supported by observations o 
experience. 

One of the results in the first sections of this paper js 
that the behavior of the ultra-long waves is sensitive ty 
even very small divergences. We may therefore e 
that if we introduce a mean divergence different from zer 
in a baroclinic model we will be able to control the ultra. 
long waves. The main problem is to relate the meap 
divergence to the parameters which carry the history of 
the flow. It means, for a two-parameter model of the 
usual type, to a height and a thickness. In the following 
simple two-parameter model we shall try to introduce such 
a net divergence in the vertical direction by relating the 
vertical mean divergence to the flow parameters in a way 
similar to the one used in the barotropic case treated 
earlier in this paper. 

The prognostic equations we are going to use will be the 
vorticity equation in the form: 


or _- Ow : 
at VY VEtN=fo dp (5.1) 


and the adiabatic equation 


5 (SY )+V- v(3 out e w=0 (5.2) 


where fy is a standard value of the Coriolis parameter. 

We shall in the following denote the 200, 400, 600, 800, 
and 1000-mb. surfaces with subscripts 0, 1, 2, 3, and 4. 
Applying equation (5.1) to the 400 and 800-mb. surfaces 
and approximating 0w/dp by finite differences we obtain 
with P=40 cb.: 

re) ¥ 
Shy, Vath= fo (, —wyo) (5.3) 

and 


oii, VEth= Jo —wy). (5.4) 


We shall in the following apply the approximate 
boundary condition o,=0, thereby disregarding mountail 
effects in the model. The practice before has been to st 
wo=0. This means that we assume that the vertical 
velocity has a zero point somewhere between the two er 
points at the surface of the earth and at the top of the 
atmosphere. Although it is true that such a zero poilt 
exists in most cases the result is that putting this zero point 
invariably at 200 mb. reduces the vertical mean divergenee 
in the layer between 200 and 1000 mb. to zero, and we Call: 
not any longer control the long waves. We shall therefore 
not put ,=0, but try to work with the layer between 100° 
mb. and 200 mb. as an “open” system. 

With »,=0 it is now convenient to add and subtract the 
equations (5.3) and (5.4). At the same time we shall 
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assume that for any quantity « we have 


a= (a +a) (5.5) 
and we shall further introduce the notation: 
a’ =5 (a,—as)- (5.6) 
Using this procedure we obtain: 
SVU tt NEV Wi =— shun (6-7) 
oF +N: VE+V’ VS +h) =a — My (5.8) 


By applying the adiabatic equation at 600 mb. the 
value for w, may now be obtained in the usual way: 


2 , i 








with 


2fo 
we am constant. 


(5.10) 





The value of w/2P may be related to the parameters 
characterizing the flow in the following way. Let us 
apply the continuity equation at 400 and 800 mb. and 
obtain 


V-Vi=— (dep). =— 5 (5.11) 
V-Vs=—(e/2p)s=F° (5.12) 
Adding these two equations we obtain: 
V . V.= 5p (5.13) 
We have on the other hand: 
ennumang 1 5P/9 d Wy 
v- V=spI,, V-Vdp=5p (5.14) 


whi Vv ‘3 the vertically averaged divergence in the 
layer between 200 and 1000 mb. According to (5.13) and 
5.14) VV. represents within the approximations of the 
model the vertically averaged divergence. Due to the 
tation (5.13) we are in the position to express wo/2P by 
‘procedure similar to the one used in the barotropic case. 

In the incorporation of the mean divergence in the two- 
pirameter model we should have in mind that we are 
ist of all trying to control the behavior of the ultra-long 
vaves. In the first approximation we are probably there- 
bore allowed to express V- V in the same way as in the 
barotropi case i.e., , 
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-%, wo=—feV > Vi=r ots (5.15) 
with 
_fit (de) aA, sie 
~a2\dp), dp mess 
With f>=10- sec.-!, o.=3 MTS units, dAldp=—"> 


and (do/dp), computed from a formula corresponding to 


(2.8) we arrive at a value of 
r=2.22 1072, (5.17) 

The assumption (5.15) means that we assume the ultra- 
long waves to have a small tilt vertically around level 
2; ie. 600 mb. It is realized that this assumption needs 
further investigation. In a paper by Eliasen [10] it is 
shown that the tilt of the long waves on the seasonal 
mean charts is not insignificant in the lower part of the 
troposphere. On the other hand, a similar Fourier 
analysis of the long waves at the levels 800, 700, 500, 300, 
and 200 mb. on a few individual days has in no case shown 
tilts nearly as great as those found on the normal charts. 
The assumption (5.15) may therefore in many cases be 
considered as a good first approximation. 

The expression (5.15) could be introduced directly in 
(5.7) and (5.8). However, diagnostic computations of 
the vertical profile of the vertical velocity described in 
section 8, show that o<o,. For simplicity we shall 
therefore neglect the second term on the right side of (5.8) 
in comparison with the first. This means only that we are 
neglecting the mean divergence in the troposphere in com- 
parison with the divergence at 400 or 800 mb., which 
represent the upper and lower parts of the troposphere. 

Granted that these approximations can be made we 
arrive at the following modified forecast equations 


ay, 
ot 


Of» 


op Te Vt N+V Ve =r (5.18) 


ar’ =O pv ; 
SEM w+ VG$) =a(Se+Ve- W ) (5.19) 


The only change which in fact has been made is that the 
term rs has been added to the first prognostic equation, 
when we compare with the equations for a model where the 
mean divergence is zero. 

It is the purpose of the following section to show that 
this small modification will control the ultra-long waves 
in the atmosphere and further to illustrate the changes 
which we may expect in the short waves using (5.18) and 
(5.19) as prognostic equations. 


6. A PERTURBATION ANALYSIS OF THE 
BAROCLINIC CASE 


Although the final proof of the applicability of a set of 
prognostic equations always is the successful forecasts 
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made from them it is of interest to make a simple analysis 
of the model in order to test it as far as the mathematical 
technique allows us to do it. The only test which can be 
made by a relatively simple mathematical technique is a 
linear perturbation investigation. The two-parameter 
models where the mean divergence vanishes have been in- 
vestigated in this way by, for instance, Eliassen [9], 
Thompson [16], and Phillips [14]. It is therefore of con- 
siderable interest to make a comparative study of the 
model described here and the corresponding model with 
vanishing mean divergence. 

The two parameters that characterize the model are the 
stream functions ¥, and y’. Let these {be {described [by 
the relations 
(6.1) 


¥o(r,y,t) = —Uzy+yre**- 


and 


v (2,y,t)=—U'y+W eto (6.2) 
where U, and U’ are constants, $, and y’ are the ampli- 
tudes, k=2z/L is the wave number, and c is the phase- 
speed (wave velocity). 

Inserting (6.1) and (6.2) into (5.18) and (5.19) one finds 
that y, and y’ are solutions to the prognostic equations 
provided 
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For any quantity a the superscript * means 


avai (6.4) 
Further 
D=[(1+q*) r* U,+(q*—r*) B*? 

—4 (1+r*)(q*—1)U”. (6.5) 


It should be mentioned that the corresponding formula 
for the usual two-parameter model with vanishing mean 
divergence simply may be obtained by putting r*=0. 
This simple relation makes a comparison between the two 
models fairly easy. 

As we are here especially interested in the ultra-long 
waves we shall first investigate what happens when 











k--0 (Lo). Considering the three terms ‘in (6.3) we 
find directly that 
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We find therefore that 


nk, if the plus sign is chosen 
lime=4 q (6.7) 
_ =“ if the minus sign is chosen 


The result (6.7) shows that both of these solutions pp. 
main finite for k->0. The interesting solution is the ong 
corresponding to the minus sign. With the numeric) 
values chosen here we find for the latter solution 


Cy-0= —B/r= —7.2 m.sec?! (6.8) 


as compared to ¢,.o=— © for r=0 corresponding to , 
vanishing mean divergence. We may therefore conclude 
that the model controls the ultra-long waves in the sense 
that the retrogression is greatly reduced. 

Next, we shall investigate in which way the baroclinic 
instability is being changed by the introduction of a mean 
divergence. Unstable solutions are found in cases wher 


D<0 (6.9 


From the expression (6.5) for D it is seen that onl) 
stable solutions are possible, when 


q*<l 
With 
g— 24X10" 12 
we find that g*<1, if Z<3.1X 10° km. Consequently all 
waves with a wavelength smaller than 3100 km. are stable 
This result is the same as for the model with vanishing 
mean divergence. 

Unstable solutions are possible if Z is greater than this 
critical wavelength. The division between the unstable 
and stable region may be found by equating /) to zero. It 
is seen immediately that a certain difference now appears 
between the models with and without mean divergence. 
In the case of no mean divergence the division bet ween the 
stable and unstable region is given by a single relation 
between the models with and without mean divergence. 
However, when we have a mean divergence the zonal wind 
U, also enters into the relation. By inspection of (6.5) i! 
is seen that the greater the velocity UV, is, the greater 1” 
has to be in order to create instability. This relation 
illustrated in figure 2, where the critical vertical wind 
shear (dU /dz). expressed now in the unit m. sec.* km” 
is plotted as a function of wavelength in units of thousands 
of km. It is seen that the curve for U/,=20 m. see. at 
higher than the one corresponding to U,=10 m. set.’ 
The curve corresponding to no zonal velocity, [/,=0 m. 
sec.-!, is also drawn as representing an extreme case. The 
dached curve in figure 2 represents (dU//dz) ¢ as : . function 
of wavelength for the model with vanishing mean diver 
gence. This curve is obtained from (6.5) by putting 
r=0. 
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Cons: lering first the ultra-long waves we find that they 
are closer to an unstable situation in the model with mean 
divergence. However, disregarding the case where 
/’,=0 as an extreme case, we find that vertical wind shears 
ys observed in the atmosphere would never make these 
waves unstable. We may therefore conclude that we have 
obtained a considerably reduced retrogression of ultra- 
long waves Without introducing any instability in this part 
of the spectrum. 

Turning our attention toward the smaller wavelength in 
the spectrum we find that the introduction of a mean 
livergence in the model changes the instability to some 
extent. We find for instance with the value of the static 
stability chosen here that the wavelength of maximum in- 
stability is around 4000 km. in both models. However, 
in the model without mean divergence a vertical shear of 
only 1.6 m. sec.-? km.~* is necessary to produce instability. 
The corresponding figures are about 3.1 and 5.8 m. sec. 
km. for 7’, equal to 10 m. sec.-' and 20 m. sec.-’, respec- 
tively. Briefly, we may therefore conclude that the way 
in which we have introduced the mean divergence into the 
model tends to stabilize the shorter waves somewhat. 

It is a general impression that a two-parameter model 
without mean divergence will develop pressure systems too 
much in cases when they actually are deepening. From 
the result of the instability analysis made here it is seen 
that a model with mean divergence will counteract this 
tendency. Actual computations will of course have to 
prove this tentative conclusion. 

It has been pointed out by other investigators (Eady 
's|, Charney [4]) that according to linear baroclinic 
instability theory the westerlies should constantly be in 
a state of instability. The result that the present model 
sems to decrease the instability makes it interesting to 
investigate whether this also is the case for this model. 
Although our instability criterion is derived for flow 
patterns without horizontal shear we shall nevertheless 
iry to compare the results with the mean state of the real 
atmosphere. The wavelength of maximum instability 
seems to be around 4000 km. We shall therefore investi- 
gate the problem mentioned above by setting Z =4000 km. 
and then computing the critical vertical wind shear as a 
function of latitude using values of 7’,; i.e., U’ at 600 mb., 
taken from a mean cross-section of the atmosphere. The 
critical wind shear may then be compared with the actual 
wind shear in the mean cross-section. For this purpose 

the mean cross-section prepared by Hess [11] has been 
used. The result of the computation is shown in figures 
stand 3b for winter and summer, respectively. A marked 
(lifference compared to earlier results is that the mean 
“mosphere is stable at almost all latitudes. The only 
exception worth mentioning seems to be north of 60° N. 
itstmimer, when the mean atmosphere seems to be a little 
unstab | 

Anot! er measure of the degree of instability which is 
present in the model is the time it takes to double the 
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FicuRE %.—Diagram giving the division between stable and un- 
stable regions for different values of the zonal velocity, Us. 
The dashed line is the corresponding curve for a model with 
vanishing mean divergence. The horizontal axis gives wave- 
length in thousands of km., the vertical axis vertical wind shear 
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FicurE 3.—The dashed curves give the averaged vertical wind 
shear in the layer between 800 and 400 mb. in the unit m. sec.” 
km.“ as a function of latitude obtained from Hess’ [11] cross- 
section. The solid curves give the critical wind shear cor- 
responding to the zonal wind as taken from the cross-section at 
600 mb. (a) Winter, (b) Summer. 
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Figure 4.—Diagram showing the vertical wind shear necessary 
to double the amplitude of a wave in 24 hours as a function of 
wavelength for different values of the zonal velocity U: The 
dashed curve gives the same for the model with vanishing mean 
divergence. 


amplitude of a certain perturbation. This time will of 
course also depend on the wind speed J, in the model 
with mean divergence. Figure 4 contains curves giving 
the vertical wind shear necessary to double the amplitude 
in 24 hours as a function of wavelength for different 
values of J, in the model with mean divergence. The 
curve corresponding to the same condition for the model 
with vanishing mean divergence is plotted for compari- 
son. This figure illustrates even more clearly than figure 
2, the decreased instability caused by the mean divergence. 
As an example taken from figure 4 it may be mentioned 
that for 7=4000 km. and U7,=20 m. sec.“ a vertical wind 
shear almost twice as large as in the model without mean 
divergence is needed to double the amplitude in one day. 

Figures 5-8 give the phase speed ¢ as a function of 
wavelength in the two models for different combinations 
of UY, and U/' The curves in figures 5a—-8a were computed 
from (6.3) with r=0, corresponding to no mean diver- 
gence, and in the figures 5b-8b with r#0. The main fea- 
ture in these figures is of course the great difference be- 
tween the curves c (corresponding to the negative sign 
in front of the square root), assuring us that no excessive 
retrogression takes place in the modified model. Figure 
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Figure 5.—(a) Phase-velocity as function of wavelength for the 
model without mean divergence, (b) the corresponding figure for 
a model with mean divergence. The dashed line in (a) is the 
Rossby speed, while the dashed curve in (b) corresponds to 
U’=0 (no vertical wind shear). Parameters: U;,=20 m. sec. *, 
dU/dz=4 m. sec. km.™. 


5 illustrates a case where the modified model is stable for 
all wavelengths, while the model without mean divergence 
actually is unstable in a band from about 3100 km. to 
about 7000 km. Note also in figure 6 that the modified 
model is unstable in a shorter interval than the non- 
modified. This is, however, not invariably so. When the 
zonal speed U, is small and dU//dz is large the modified 
model will be unstable in a broader band of wavelength. 
This fact is illustrated in figures 7a and 7b, corresponding 
to 17,=10 m. sec.” and dU /dz=8 m. sec. km.-'. Figures 
8a and 8b finally illustrate a case where U7,=10 m. see. 
and dU /dz=4 m. sec. km.. In these diagrams the band 
of unstable waves is about the same. 

Figures 5a-8a contain for comparison the speed of 
waves in a non-divergent, barotropic atmosphere, the 
Rossby speed, while figures 5b-8b contain the wave speed 
corresponding to dU/dz=0. For large values of the 
wavelength the speed for d77/dz 0 and dU’/dz=0 are 
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FieurE 6.—Same as figure 5, but with parameters: U,—20 sec.™’, 


dU/dz=8 m. sec. km.~. 


about the same, while the speed of the waves in a baro- 
clinic atmosphere is somewhat larger than in an atmos- 
phere with no vertical wind shear. 


7. ON THE VERTICAL VARIATION OF STATIC 
STABILITY 


The evaluation of divergence in sections 2 and 5 depends 
toa large extent on the vertical variation of static stabil- 
ity. We have in the earlier sections used a variation 
described by (2.8), where the static stability varies 
inversely to the square of pressure. 

The tropospheric part of a mean atmosphere is usually 
characterized by a certain lapse-rate, y= —27'/2z. The 
most straightforward way to compute the static stability 
=~a°ln@/Op as it appears in models for numerical inte- 
gration would therefore be to assume a certain value for y 
id then compute the variation of o. It is, however, of 
great convenience in many problems to have o as a rather 
simple function of pressure. In this section we shall show 
that a formula of the form (2.8) will describe the vertical 
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Same as figure 5, but with parameters: U,—10 m. sec.”, 
dU/dz=8 m. sec.” km.”, 


FIGURE 7. 


variation of temperature in the atmosphere to a good 
approximation. Let us write (2.8) in the form: 
a = 
fame. 4 (7.1) 
Py 
where @ is a constant determined in such a way that o at 
some pressure level, for instance 600 mb., is equal to a 
standard value determined from a mean atmosphere, Let 
us next try to find 7’=7'(p) corresponding to (7.1) by 
writing 


dIne R/RT aT ‘ia 
On o_ 


oc=— —{ (4.2) 


op pe,p dp 


we find that the temperature has to satisfy the following 
differential equation of the Euler type 


dT R a “ 
oh ct —— 7 
p ws, T+7 0. 7.3) 
The solution to (7.3) is 
R 
T= T(p)= Frito + Cp. (7.4) 
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FieurE 8.—Same as figure 5, but with parameters: U.=10m. sec.", 
dU/dz=4 m. see.” km.. 


The arbitrary constant C can be determined by fixing 
the temperature at a certain pressure, say, 600 mb. We 
obtain then 


L. 

7=7{ P)= Fr rt( Tif o> )( 2) (7.5) 
The temperature distribution given by (7.5) is plotted 
in figure 9 as the full curve. The dashed curve gives the 
temperature distribution in the standard atmosphere, 
while the dotted curve is taken from Defant and Taba [7] 
corresponding to the air mass between the polar-front jet 
and the subtropical jet. 
follow each other quite closely, and we may therefore con- 
clude that the assumption (2.8) represents mean conditions 

in the troposphere in the westerlies quite well. 
For use in the next section we shall need the tempera- 
ture distribution in the stratosphere. 


It is seen that the three curves 


For simplicity we 
shall assume that the stratosphere is isothermal with a 
temperature corresponding to the one obtained from (7.5) 
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Figure 9.—Temperature as a function of pressure corresponding to 
a variation of static stability inversely proportional to the 
square of pressure (solid curve). The dashed curve gives tem- 
perature as function of pressure in the standard atmosphere, 
while the dotted curve is the temperature as function of pres- 
sure in the air mass between the polar jet and the subtropical) 
jet as given by Defant and Taba [7]. 
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Figure 10.—Static stability as a function of pressure in tropo- 
sphere and stratosphere, when the stability is inversely propor- 
tional to the square of pressure. 





for p=p.=200 mb. The stability will then vary in the 
following way: 


RT; 1 4 
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(7.6) 
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The stability profile as determined from (7.1) and (7.6) 
is shown in figure 10. It is seen that the static stability 
varies greatly in troposphere and stratosphere and further 
that the tropopause is characterized by a jump in st ability. 
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3, SOME EXAMPLES OF VERTICAL PROFILES OF 
THE VERTICAL VELOCITY 


In section 5 we neglected the vertical velocity at 200 mb. 
in comparison with that at 600 mb. As the 200-mb. sur- 
face usually in the westerlies is in the stratosphere one 
can from many considerations justify this neglect. In 
order actually to measure the ratio between the two 
quantities certain diagnostic computations of the vertical 
profile of » have been performed. These computations 
were made with a model slightly different from the one 
described in this paper, because special attention was paid 
to the influence of the stratosphere on the vertical profile 
of. ‘The w-equation was obtained from the vorticity and 
adiabatic equations in the forms 


or ] or Ow 
sit VEtA +o = Cth dp (8.1) 
and 
od (dv wae. 
ey] Se )tV Vv (S¥)+50=0 (8.2) 


by the usual procedure. In this way we arrive at an 
equation for the vertical velocity in the form 


" o* ; 
(Vy+f) ap tere p dp OnE sy?) (8.3) 
with 
’ i . ‘eaK 2 . oy ¥ 
M@wP)=95 (V°Vn)—V (Vv Vv 3) (8.4) 


In the troposphere we assume that 


V=V.+A(p)V’, c= 7 (8.5) 

and in the stratosphere that 
b 
V= B(p)WVo, = 7 (8.6) 


The functions A(p) and B(p) were in fact determined 
insuch a way that they are consistent with the assumption 
that Vo—0 (see Berkofsky [1]) and further such that the 
wind is continuous at p=po=200 mb. (see fig. 11). When 
the approximations (8.5) and (8.6) are introduced in (8.3) 
we arrive at two equations, one for the troposphere and 
one for the stratosphere. 

As we are here interested in the vertical variation of , 
we shall replace Y*%o by — fo, where £ is an inverse 
ieasure of horizontal scale. By doing this we restrict 
ourselves to a consideration of essentially one scale. 8 was 
here chosen as corresponding to the wavelength of maxi- 
mum instability (Z=4000 km.). With this restriction 
the two equations obtained from (8.3) reduce to ordinary 
differential equations in », which then should be solved 
using ajypropriate boundary conditions. 
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Figure 11.—The functions A(p) and B(p) characterizing the 
variation of the horizontal wind with pressure in the model 
used for diagnostic computations of vertical velocity. The 
dashed curve in the lower figure corresponds to the case where 
A(p) is a linear function of pressure. 


We shall use o=0 for p=0 and p=p,=1000 mb. 
Furthermore, we need an internal boundary condition for 
p= po=200 mb., where the solution to the stratospheric 
equation should be matched together with the solution to 
the tropospheric equation. The level p=p,=200 mb. 
represents in the present considerations the tropopause as 
well as the level of maximum wind. These two levels 
need not in all cases be the same in reality, but synoptic 
investigations (Defant and Taba [7]) show that they are 
close to each other, usually with the level of maximum 
wind a little lower than the tropopause. We way, 
therefore, at the level p=p, apply a boundary condition 
natural for a discontinuity surface for stability. Such 
boundary conditions have been treated by Lowell [12]. 
The vertical velocity should be continuous at the tropo- 
pause and further should (d/dp)s~(0w/dp)r which 
means that the first derivative in the p—direction should 
be continuous. These two conditions have therefore been 
applied at p= po. 

The solution to the two equations obtained from (8.3), 
with the boundary conditions mentioned above, is quite 
laborious. For the general case it has not been possible 
to find analytical solutions. A numerical integration of 
the two equations has therefore been performed. The 
method used in the numerical solutions has been the 
Runge-Kutta method, which is most convenient in cases 
of internal boundary conditions. Examples of solutions 
are presented in figures 12 and 13. In these cases the 
stream functions at 600 mb. and for the thermal flow 
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Figure 12.—The vertical velocity as a function of pressure in dif- 
ferent points of a baroclinic wave, where the thermal stream 
function is lagging % of a wavelength behind the stream func- 
tion at 600 mb. (see upper part of the figure). The dashed curve 
in the upper part of the figure gives the vertical velocity through 
the baroclinic wave at 600 mb. 
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Figure 13.—Same as figure 12, but with phase-difference of 4 of 
a wavelength between the thermal stream function and the 
stream function at 600 mb. 


between 800 and 400 mb. were specified as simple sinus- 
oidal waves having a certain phase difference. Figure 
12 corresponds to the case where the thermal field is lag- 
ging 14 of the wavelength behind the field at 600 mb., 
while figure 13 has a phase difference of 14 of a wave- 
length. All constants were computed for eight different 
points. The computation gave then the eight vertical 
profiles represented in the figures, which in the upper part 
contain a schematic picture of the temperature and pres- 
sure wave and also the horizontal profile through the 
baroclinic wave of the vertical velocity at 600 mb. On 
each of the profiles is given the level of maximum vertical 
velocity and the levels where » has a zero point different 
from those at the boundaries. It should further be men- 
tioned that each profile is normalized in such a way that 
the maximum value is unity. 
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Of importance for the discussion in the preceding gee. 
tions are the following features: 

(a) A zero point for the vertical velocity appears some. 
what higher than the level of maximum horizontal! wing, 
The pressures of the zero points vary systematically 
through the wave, and we may state that the values of 
» obtained above the zero point are extremely small, 

(b) wo at 200 mb. is not exactly zero, but it is true that 
© Kw, Where o, is at 600 mb. This characteristic justifies 
the approximation made in section 5. 

(c) The non-divergent level varies systematically 
through the baroclinic wave with the result that this level 
does not coincide with a pressure level. 

(d) It isa limitation in the present calculation that the 
tropopause level has been assumed at a constant pressure 
(200 mb.). In reality we may therefore expect a larger 
variation of the level of nondivergence. 


9. SUMMARY AND CONCLUSIONS 


Sections 2-4 contain a derivation and discussion of the 
prognostic equation for the 500-mb. flow where attention 
is focused on the term which controls the ultra-long waves, 
It is shown that the numerical value of the coeflicient is 
related to the vertical variation of the horizontal wind and 
the corresponding variation of the static stability. The 
numerical value determined in this way agrees well with 
the one obtained empirically by Cressman [6] and Bolin 
[2]. 

Section 5 generalizes the results obtained in sections 
2-4 to a simple two-parameter model. It is shown that 
the simplest way to control the ultra-long waves in a baro- 
clinic model is to introduce a mean divergence in the 
troposphere. The mean divergence may in the first ap- 
proximation be estimated barotropically corresponding to 
the well-known fact that the mean motion in the atmos- 
phere is almost barotropic. 

The influence of an existing mean divergence in the 
atmospheric layer under consideration on the baroclinic 
instability is investigated in section 6. The main result 
is that the mean divergence stabilizes the shorter waves 
to some extent, beside the effect of decreasing the retro- 
gression of the ultra-long waves considerably. An appli- 
‘ation of the stability criteria derived in this section to 
a mean cross-section for the atmosphere shows that the 
mean atmosphere is baroclinically stable at almost all 
latitudes. 

A justification of certain approximate relations used in 
the earlier sections is given in sections 7 and 8 of whieh 
the former is concerned with the vertical variation of the 
static stability, while the latter gives certain diagnosti¢ 
results of computations of vertical profiles of the vertical 
velocity. 

Any model designed for numerical prediction of trope 
spheric flow pattern must contain an effect which at least 
expresses the quasi-stationary behavior of the ulira-long 
waves. According to results obtained by Burger [3] We 
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cannot expect to express more than this as long as we use 
the vorticity equation as the prognostic equation. One 
jossible way of controlling the ulfra-long waves is demon- 
strated in this paper, where we have made use of the 
observed barotropic character of these waves. Certain 
modifications of the behavior of shorter waves is also 
obtained by this procedure. It is quite likely that a more 
sophisticated introduction of the mean divergence, which 
seems to be the important quantity, into a tropospheric 
model inay change the behavior of the Rossby type of 
waves still more. 
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CORRESPONDENCE 


Comments on “The Hurricane Season of 1958 


aa 


EUGENE W. HOOVER 


District Forecast Office, U.S. Weather Bureau, Washington, D.C. 
May 19 1959 


In a recent article [1] the Staff, Weather Bureau Office, 
Miami, Fla., examined the interaction between hurricanes 
Helene and Ilsa, 1958, and suggested that their relative 
motion had been counterclockwise. Reexamination of the 
data shows that the relative motion was clockwise and that 
sume changes should be made in table 1 and figure 7 of the 
article to take into account the following items: 

1. The motion of Ilsa relative to the midpoint must be 
exactly equal in magnitude but opposite in direction to the 
motion of Helene relative to the midpoint. This follows 
because by definition the midpoint remains the same dis- 
tance from each of the storms. Thus, in table 1 in the 
columns headed “Ilsa-- Midpoint” and “Helene—Mid- 
point,” the two courses at each time should differ by 180° 
and the two speeds should be equal. 

2. The speed of the midpoint cannot exceed the arith- 
metic mean of the speeds of Helene and Ilsa. Thus, in 
table 1, under the column headed “midpoint,” the speed of 
9.6 knots at 0000 emr on the 27th is not consistent with 
Ilsa’s motion of 335°, 8 knots and Helene’s motion of 305°, 
10 knots. 


The midpoint, whose motion is the vector mean 


of the motions of Helene and Ilsa, should have a speed less 
than 9 knots. 

3. The relative motions of Helene and Ilsa should not be 
considered with reference to a west-east line, as in figure7 
of the article, but with reference to the line connecting 
their centers at the beginning of the period, i.e., a line 
oriented approximately northwest-southeast. The cor- 
rection of figure 7 by properly reorienting the dashed lines 
shows that the centers of the two storms had a clockwise, 
converging relative motion during the last 36-hour period 
covered by table 1. 

It is not suggested here that the concept of interaction 
between vortices, which calls for counterclockwise relative 
motion, is not valid, but in this particular case, during the 
36-hour period, other factors were exerting a greater in- 
fluence to produce a net clockwise relative movement of 
the hurricanes. 
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NOTE ON FREQUENCY OF HIGH WINDS OVER THE UNITED STATES 


H. E. LANDSBERG AND B. RATNER 
Office of Climatology, U.S. Weather Bureau, Washington, D.C. 


[Manuscript received April 28, 1959; revised May 20, 1959] 


Problems of jet aircraft operations and radioactive fall- 


out have created much interest in high-speed upper-air 
currents. Also, ever since Rossby’s [4] classical paper 


there has been a lively interest in the dynamics of the jet 
gream and associated models of atmospheric circulation. 

Some earlier climatological assessment of the mean posi- 
tion and strength of the jet stream is contained in the 
monograph on the jet stream by Riehl et al. [3]. In recent 
years more reliable observations of high-speed winds have 
ecome available. Since 1957 the National Weather 
Records Center has kept special records on punch cards of 
ull very high wind speeds at United States stations. 

For purposes of quick description we set an admittedly 
wbitrary limit of 100 miles per hour and had a machine 


compilation made of the frequency of observations equal 
to or above this value. This collection covers 36 stations 
with the best upper-wind-observing equipment (either 
GMD-1 orGMD-1A). It is unlikely that many observa- 
tions there were missed because of limiting angles. 
Although only the data for 1957 and 1958 were available 
for summarization, the results are sufficiently interesting 
to warrant a preliminary note. The data are presented 
in the form of three charts. Lines of equal frequency of 
occurrence indicate the distribution of high-speed winds 
in January, July, and for all observations in the two years. 
Figure 1 shows the midwinter conditions. The axis of 
highest frequencies coincides well with the corresponding 
average position of the jet stream given by Namias and 
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Figure 1.—Percentage of soundings with winds = 100 m.p.h. January 1957-58. 
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Clapp [2]. It is concentrated just north of the Gulf of 
\fexico inland and extends eastward to Cape Hatteras. 
jt that point on the Atlantic coast not less than 77 per- 
ent of the ascents showed winds 100 m.p.h. in the two 
Januarys. 

During July, shown in figure 2, the maximum of high 
yind frequencies occurs Just north of the Great Lakes and 
extends eastward to the St. Lawrence Valley. A little 
over 30 percent of the observations in this area show high 
windsaloft. The core seems to be a little farther east than 
in the corresponding average jet stream position of 
Yamias and Clapp [2]. This difference may well be 
used by the very short record in our series, while their 
jata were based on calculations from mean upper-air 
charts. Notable in figure 2 is the almost complete absence 
of winds above 100 m.p.h. south of 35° N, latitude. 

The mean annual pattern shown in figure 3, although a 
mixture of all seasons, has the virtue of having the highest 
number of observations. These average about 730 per 
station in the chosen network. The meridional pattern is 
quite pronounced. There is no zonal high-speed band, at 
least for the frequencies of winds 100 m.p.h. These high 
speeds are most common over the northeastern area where 
they reach over 40 percent of all observations. This maxi- 
mum is located at about 42° N. latitude. 

Interesting also is the narrow zone of increased high- 
seed wind frequency just east of the Rocky Mountains, 
followed immediately to the east by a minimum over the 
Great Plains. Also notable is the low frequency of high- 
speed winds over the central California coast. An effect 
of the High Sierras is possible but can not be documented 
by the present network. 
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The distribution east of the Rocky Mountains can be 
explained by the presence of a dynamic trough in the 
lee of the mountains. Boffi [1] discussed this in connection 
with his analysis of the stream field produced by the 
Andes. The marked asymmetry of west and east coasts 
has already been commented upon by Riehl et al. [3]. 
The northeast maximum in the annual pattern as well as 
the winter maximum centered over the Carolinas and 
Georgia would seem to require injection of additional 
energy into the basic dynamic system envisaged by Rossby 
[4]. Previous writers have belittled the idea of thermal 
contrasts being associated with the development of wind 
patterns. Perhaps in conjunction with conversion of 
latent heat into kinetic energy the thermal contrasts are 
likely to contribute both to the frequency and the specific 
latitudinal distribution of the high speed jet streams. 
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A GRAPHICAL METHOD FOR COMPUTING HORIZONTAL TRAJECTORIES 
IN THE ATMOSPHERE* 


HUGO V. GOODYEAR 


Hydrologic Services Division, U.S. Weather Bureau, Washington, D.C. 


[Manuscript received November 19, 1958; revised April 30, 1959] 


ABSTRACT 


Graphical solutions of the horizontal equations of motion are developed for computation of dynamic tra- 


jectories. 
is suggested, and graphical application is outlined. 


1. INTRODUCTION 


The purpose of this paper is (1) to study the quantita- 
tive relationships between the changes of the pressure field 
encountered by a parcel and the changes in its velocity, 
and (2) to apply the relationships to the construction of 
horizontal trajectories. 

Trajectories of air parcels have been computed by vari- 
ous investigators using the actual winds, geostrophic 
winds, gradient winds, constant-absolute-vorticity con- 
siderations, and combinations of the above. One of the 
first contributions to a dynamic approach was that of 
Machta [1]. Recently additional work on “dynamic” 
computation of trajectories has been published by Taka- 
hashi et al. [2], Franceschini and Freeman [3], Wobus 
[4], the Air Weather Service [5], and Hubert [6]. The 
methods to be presented in this paper give exactly the 
same results obtained by Wobus’ method, and by the Air 
Weather Service method for the frictionless case. Fran- 
ceschini and Freeman’s method differs slightly. Hubert’s 
method also differs, his approach being based entirely on 
hourly forecasts of geopotential or stream function fields 
by numerical methods. The handling of the frictional or 
damping terms herein differs considerably from that of 
Wobus. The other investigators have applied the methods 
only in the higher troposphere and the stratosphere, while 
in this paper the methods are applied near the ground. 


2. THEORY 


The problem is to find the trajectory of a parcel of air 
in the frictionless atmosphere when its initial velocity is 
known (or can be closely approximated) and the contour 
or pressure field is known or can be approximated during 
a later interval of time. The simplified equations of 
motion in the atmosphere without friction are 


*This paper was written as part of the work under P.L. 71, 81st Congress 
under funds transferred from the Corps of Engineers. 


Expedient techniques for using these solutions are discussed. 


An approximation for the frictional effect 


du 


a 








fog=fo' (1) 


and 


dv , 
Tp It Sle — fu (2) 


where u and v are the east-west and north-south compo- 
nents of the wind, f is the Coriolis parameter, the subseript 
g denotes the geostrophic wind, and wu’ and v’ denote the 
geostrophic deviations; i.e., u’ =u—u, and v’ =v—»,. 

Haurwitz [7] shows these equations have the simple 
geometric interpretation that the wind is accelerated along 
a line at right angles to the geostrophic deviation vector 
V’ and in a right-hand or clockwise sense, as illustrated in 
Note that a must always be directed to the 
right of V’ in the Northern Hemisphere. 

For computing the change in the wind during a certain 
interval of time, one may use the average geostrophic wind 
on the parcel during that interval, thus assuming that itis 
constant in magnitude and direction. Then, it follows 
from equations (1) and (2) that 


figure 1. 


du’ 

— )/ (3 
Ti fi ) 
dy’ 

—— 4 
i fu (4) 


during the interval. Now, from figure 1 

u’=V’ cos 6’ (5) 
v’=V’ sin 6’ (6) 
where V’=|VW'|=vu'2+0". Use of (5) and (6) to trans- 


form equations (3) and (4) into a pair for V’ and # gives 


>. 
dt — 





0 (7) 
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Ficure 1—Geometric relationship between the acceleration vector 
d V/dt and the geostrophic deviation vector V’. 
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These equations mean that the ageostrophic wind vector 


is constant in magnitude and rotates anticyclonically with 
aperiod of one-half a pendulum day. 

Thus, during the time interval in which the geostrophic 
wind is considered a constant (equal to its average), the 
change in the wind vector can be represented by the rota- 
tion of the ageostrophic wind vector, keeping constant 
length, clockwise at a rate equal to the Coriolis parameter. 

Now, letting Vi=initial wind, and V;=final wind after 
time interval + 


Vi=Vit fr’ Xk (9) 


where the k is the vertical unit vector. 


3. GRAPHICAL SOLUTION 
THE ACCELERATION TERM 
ne (he vector V’ is constant in magnitude, and since 


@ 7, we may devise a simple graphical solution to 


evaluate the acceleration term. 
In figure 2, Vi, V,, and V;’ represent the initial wind, the 
geostroplic wind, and the initial geostrophic deviation. 


MONTHLY WEATHER REVIEW 189 


fe ; 














Oo 


Ficure 2.—Geometric relationships among the initial wind V;, the 
geostrophic wind V,, the initial geostrophic deviation V;’, the 
final geostrophic deviation V;’, and the final wind V;. 


Then V’ rotates anticyclonically at the rate f. At the end 
of a time increment +r, a simple rotation with an angle 
6=fr gives the final geostrophic deviation V;’. Then a 
vector from O to the end point P of V;’ gives the final 
wind V;. 

In the application of the relations shown in figure 2 to 
find the trajectory of a parcel of air during the next hour, 
the following steps are taken : 

(1) Beginning at point O lay off on the weather map 
a distance OS in nautical miles equal to the observed wind 
V; in knots. 

(2) Measure the average geostrophic wind V, along 
the path OS. 

(3) Lay off OR in nautical miles equal to V, in knots. 

(4) With pivot of compass on terminal point R strike 
an are SP of angular magnitude fr clockwise from ter- 
minal point S. 

(5) The distance OP in nautical miles equals the final 
velocity V; in knots at the end of one hour. 

(6) The trajectory during this first hour is along a vector 
which is the mean between V, and V,. A simple means 
of obtaining this average will be treated later. 

The striking of the are in step (4) explains use of the 
term “‘are-strike”’ as a name for the method. 


NOMOGRAM TO AID IN DRAWING TRAJECTORY 


Trajectories could be computed as described above by 
aid of compasses and a straight-edge, but since this in- 
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volves measuring distances and angles either on the ana- 
lyzed chart or on a working chart it is a rather slow process. 
The work can be speeded up by the use of nomograms. 
In figure 2 it is obvious that because the angle of strike 
is fixed for any given latitude and time increment and 
because V,/\=!V,'! point P will lie along a line a through 
point S whose angle 8 with line m is dependent only on the 
angle of strike, 6. Specifically, since triangle RSP is 
isosceles 
p—180°—0 sai 
To find the position of P along line a, draw a line n 
through P parallel to m and construct a perpendicular 
line k from R to line n. Draw line SN. Line SN helps 
to locate point P because it can be demonstrated that its 
angle y is also dependent only on the angle of strike; 
specifically an examination of figure 2 shows 


y=arctan (sin 6) (11) 


It is possible then to find point P by knowing merely 
the positions of the end-points of V, and V,, and the lati- 
tude and time interval. 

The next item to consider is the average velocity during 
the time interval involved in drawing the trajectory. The 
velocity with which the parcel is carried forward should be 
an integrated average V, from the initial velocity V, to 
the final velocity V;. 

If V, is not too different from V,, 

] 
Vix5 (VitV~) (12) 

If V and V;, are very different, this approximation 
becomes less accurate and then it is necessary to approach 
this from a more rigorous standpoint. Wobus [4] does 
this averaging by means of a nomogram called a Trajectory 
Computer. Experience dictates that the shorter the time 
interval used the more accurate the trajectory will be 
defined—this is particularly so where the contour or 
mressure fields are changing rather rapidly in time or space. 
We must look for a system to construct our nomograms 
rapidly forsany time interval and latitude desired. 

In figure 3 is reconstructed the portion of figure 2 neces- 
sary to aid in finding the average velocity V,. S and R, 
which are the terminal points of V; and V,, define V’ at 
time t=0. At some later time ¢t, V’ has rotated through 
angle @ and point P is then the terminal point of V. 

It may be seen from figure 3 (see also equations (5) and 
(6), noting 6=2x—8’) that 


u*—Vi—V, cos 0 (13) 
v*¥—V{ sin 6 (14) 


p* 


where u* and 
V* from § to P. 

From equation (8) and the relation 6=2—8@’, it follows 
that @6=/t; thus equations (13) and (14) can be integrated 


are components of the velocity vector 
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Figure 3.—Geometric determination of terminal point T of tra. 
jectory from vector relationships shown in figure 2 


for the time interval 0 to r, giving 


yr? 


t;= Vi'r—- . sin @ (15) 
f 
WF (1—cos 6) (16) 
Since V,—V,=V*, equations (15) and (16) give the 


coordinates of the parcel relative to point S after a time 
interval 7. If SA is the line through S and the point 
(r;, yy), making the angle @ with line m, then from figure 
3 and (15) and (16): 

Yr_ __|—cos 6 


: (17) 
ry 6—siné 


tan a= 
Now since the average velocity WV, for the particular 
arc-strike, @, has the components z,/r and y;,/r, it follows 
from (17) that the end point of the vector V, lies along a 
line which makes angle @ with line m, and @ is dependent 
only on the angle of strike; specifically, from (17): 


1—cos 6 (18) 


a=arctan 
é—sin @ 


Similarly, by integration of equations (5) and (6), it 
can be shown that z,/r and y,/7 lie on a line through 
point R that makes angle 6/2 with line m. Thus in figure 
3, ¥=0/2. 

With equations (10), (11), and (18) it is possible to 
construct quickly a nomogram to compute trajectories a 
any latitude and for any time period. Example: to com 
struct a computing nomogram for the 30°N. latitude and 





for intervals of 2 hours. The angle of strike @=/t-—-42sin¢ 
=30° where 2 S00 » the angular speed of the earth’s10- 

: 180°—@ 180° —30° _weo 
tation. From equation (10), 8=——3——-= > -=15. 


From equation (11), y=aretan sin 30°=26.6° 
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FIGURE 4. Nomographic determination of terminal point T of trajectory. The lines emanating from point S at angles a, 8, and 
are fixed for constant latitude and time interval. 





From equation (18 —arct 1—cos 6 
jua (18), a=arctan a 
—_ , ° 
agian 2 ~cos BV =80°. 


sin 30° 


On & }ece of graph paper of suitable size draw three 
lines wit! angles of 75°, 26.6°, and 80° all emanating from 


the origin S as in figure 4. For convenience the trajectory 
is best drawn on transparent paper. Lay this over the 
map and mark reference points so that the same relative 
position of paper and chart may be easily found. Mark 
the point where the trajectory is to be started. (The 
initial wind here should be known or closely approxi- 
mated.) From this point lay off a line twice the vector 
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Figure 5.—Example of construction of trajectory for model pres- 
sure field with observed wind V,, mean geostrophic wind V,, angle 
of arc-strike 6, and final wind V;. All wind vectors are drawn 
to lengths corresponding to 3-hr. displacements so that the 3-hr. 
trajectory is given by the mean wind vector Va. 
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Ficure 6.—Overlay for rapidly constructing trajectories at latitude 
85° for a polar stereographic map base of scale 1:10". 
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Vi, since the time interval is two hours. 

As a first approximation assume this is the act val path 
taken by the parcel and measure the average geostrophic 
wind along this path. References [3], [4], and | 5] haye 
some good discussions on several methods of ascertaining 
the average geostrophic wind and the relative merits of 
these methods. 

From the initial point of V; lay off twice V,. Now lay 
the transparent work-sheet on the computer so that the 
terminal end, S, of V; falls on the origin and the terminal 
end, R, of V, falls along the abscissa. This is shown jn 
figure 4. From R follow the ordinates upward to line y, 
then along the abscissa to line 8 and mark point P. A Jine 
from O to P is the final velocity V; at the end of 2 hours 
With a straight-edge along R and P mark point A on line 
a and mark a point T halfway from A to S. 

The line OT now represents the average velocity V, 
during the first 2 hours, and since it is automatically 
doubled it is also the trajectory of the parcel for the first 
2 hours. V; now becomes the initial velocity for the next 
2-hour period and the process is repeated. 

If the pressure or height gradients are not nearly wni- 
form along and near the path along which the geostrophiec 
average was measured, second and perhaps (rarely) third 
approximations should be made. 

If a second approximation is necessary in the example 
of the computation of the trajectory for the first 2 hours, 
place the work-sheet back on the chart and measure the 
geostrophic wind average along V.. Lay off the new aver- 
age V, and re-compute V; and V.. Do not change Vi, 
however. 


RAPID APPROXIMATE CONSTRUCTION OF TRAJECTORIES’ 


A method for the rapid construction of frictionless tra- 
jectories, based upon equation (12) is presented below. 
Although the method applies to any latitude or map base, 
an example is given for 35° N. latitude and a map base of 
1:10", polar stereographic projection. Figure 5 is a model 
pressure pattern showing isobars at a 1-mb. interval. A 
trajectory is to be constructed at 3-hour intervals begin- 
ning at point O. In order to remove partly the effects of 
surface friction, the so-called “gradient” wind is used 4 
the initial wind at point O. A sheet of transparent paper 
is placed on top of the map and this wind is projected for 
3 hours, or in this case 90 nautical miles, downstream to 
form the vector Vi. 

Next, the mean geostrophic wind direction for the 
3-hour period is determined at the midpoint of V; (poi! 


=). This same direction is laid off from O to form the 


line g. The 3-hour mean geostrophic wind speed is meas- 
ured using a transparent overlay illustrated in figure 6 
The upper part of the overlay gives the distance moved by 
a parcel moving with the mean geostrophie wid. The 
scale on the left is in nautical miles. The straight and 
curved lines are spaced at intervals of 20 nautic! miles. 


*The method was developed by Mr. K. Peterson. 
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In ord-r to obtain the 3-hour mean geostrophic wind 
speed, | ‘ace the scale on the map so that the straight lines 
on the overlay are parallel to the line g and so that point 


Vi 
9 





"Ties on the base line. The point——should be approxi- 


mately mid-way between the point where an isobar touches 
the edge of the overlay at zero and the point on the base 
line 2 mb. to the right. Holding a pencil over the point 
which is 2 mb. to the right on the baseline, move the scale 
so that this point is coincident with point O and the 
straight lines on the scale are parallel to line g. Place a 
dot at the edge of the scale on line g. This determines the 
ector V,, the distance a parcel starting at point O will 
travel in 3 hours if the wind is geostrophic. If a 6, 2, or 
{14-hour travel distance is desired, use the overlay to 
measure across 1, 3, or 4 mb., respectively. Such an over- 
lay can be prepared for any latitude, map base, map pro- 
jection, and isobar interval with the aid of a standard 
geostrophic wind scale and a distance scale. If overlays 
are prepared for every 5° latitude, the error at latitudes 
between overlays is negligible. 

Place the bottom point of the overlay at point V, and let 
point V; touch the left edge of the scale. Using the ares 
as a guide, go to the right (anticyclonically) and place a 
point at the right edge of the scale, locating point V;. The 
vector V,; represents the actual wind after 3 hours. Bisect 
the chord from V; to V; to obtain the point V,. (It can be 
shown that V, should be between the midpoint of the are 
and the chord, but the error is negligible). The vector V. 
is the mean actual wind given by equation (12) for the 
time period; it represents the trajectory for the first 3 
hours. 

The angle at the bottom of the overlay is obtained from 
the relation 6=fr; here fr=30°(sin @)r where ¢ is the 
latitude and + is the are-strike interval. At @¢=35° and 
with a 35-hour are-strike, 

6=30° (0.574) (3) =51.8° 
The above equation can be used to obtain arc-strike angles 
for any latitude and are-strike interval. 

In order to obtain the next 3-hour trajectory, transfer 
the transparent paper to the next map, in this case, 3 hours 
later. Take the vector V, and translate it so that the 
initial point O is coincident with the point V.. This new 
vector is used as the actual wind for the next 3 hours. The 
remainder of the technique is the same as before—draw 
the mean geostrophic wind direction, use the overlay to 
determine the mean geostrophic movement, and strike the 
are to obtain the new 3-hour final wind and mean wind. 
This procedure can be repeated indefinitely, but since the 
trajectories are frictionless, and since their accuracy is de- 
pendent on the accuracy with which the geostrophic wind 
field, both present and future, can be delineated, they must 
be use! with care when projected for much more than 12 
hours. While the large-scale features might be forecast 
with si ‘ficient accuracy, the fine-grain structure might not, 
and an result dependent on the latter will suffer. 
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4. FRICTION EFFECT 
EFFECT OF FRICTIONAL ACCELERATION ON TRAJECTORIES 


Without friction, once a trajectory is started it may be 
carried on step by step, chain fashion, for as long a period 
as the contour field is known or can be approximated. 
However, experience has shown that this tends to give 
unrealistic results as inertial oscillations dominate the 
motion. 

The introduction of some damping effect will give more 
reasonable trajectories but will not, of course, eliminate all 
the difficulties. The major effect of the viscous or fric- 
tional forces is to slow down the parcel speeds, and this 
neglect leads to systematic errors. For this reason, con- 
siderable effort has been directed to devising a method of 
applying frictional corrections. The simplest hypothesis 
regarding the frictional forces is that they are directed 
against the wind in direction and are proportional to the 
first or some higher power of the speed. It is suggested 
that a closer approximation to reality can be obtained if 
a correction is made using such a relationship with the 
constants of proportionality determined statistically from 
a comparison between trajectories computed by dynamic 
and kinematic methods for a large number of observations. 
The suggested procedures are described here; evaluation 
of the empirical constants will be reported at a later date. 
Doubtless, the more modern theories of the frictional 
forces can be incorporated by more elaborate procedures. 


THEORY 


If the acceleration due to friction of an unspecified 
nature is defined such that it is proportional to a constant 
power of the speed, the equations of motion may be 
written. 


. er . 
an —ku 
om —fu’ —kv" 


These equations can be solved graphically if n is known. 
probably n is a function of the speed itself, but for sim- 
plicity a value of 2 was assumed. 

The frictional effect, as defined above, was added at the 
end of each time interval to the computations of the 
trajectory by the methods already described. This in- 
volves finding the time-average of the friction term kV? 
during the interval of time 7 in which V, changes to 
Vy, where V; is the final speed after the friction effect has 
been added. 


oe f, "EV 'dt=V,—Fr (19) 


where F is average frictional acceleration, Vr is speed at 
time r. The average acceleration a during the interval r is 
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FRICTION NOMOGRAM 
k =.003 +t=1 hour 
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7.—Friction nomogram based on equation (22) with 


k=0.003, r =1 hr. 


FIGURE 


Putting V=V,+@t and using equation (20) to eliminate @, 
we can evaluate the integral in the middle member of (19). 
This gives 


T loa ae 
[irae [rn 2s] an 
0 e 
Putting (21) into (19) gives 
Vr=V,—kr [Vit ES a (22) 


GRAPHICAL SOLUTION OF EQUATIONS 


Equation (22) can be solved graphically. See figure 7. 

In equation (22) kr occurs outside the brackets in the 
last term. This suggests that once a system of isolines 
of friction deceleration values is calculated this same chart 
can be used for other values of % and + by the simple 
expedient of multiplying the values of the isolines by an 
appropriate factor. This expedites the testing of different 
values of & in different situations, either actual or models. 
To use the friction nomogram, figure 7, enter the abscissa 
with the value of the initial wind speed V; and the ordi- 
nate with the value of the final wind speed V; obtained 
by the are-strike technique. At the point thus obtained 
interpolate between curved isolines. This interpolated 


value is subtracted from V; to yield V»y, the final speed 
after reduction due to friction. 
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APPLICATION TO TRAJECTORIES 


As previously mentioned, the friction effect is added 
stepwise at the end of each trajectory computation. This 
can be considered only an approximation. Nevertheless, 
if the method proposed is to be useful, it must be able to 
predict certain singularities in the field of motion. For 
example, any wind, initially not in balance, in a straight 
uniform pressure gradient field will eventually approach 
a value such that the accelerations on any parcel of air 
total zero. At this point there is a balance among the 
forces due to the Coriolis effect, the pressure gradient, and 
friction. 

Furthermore, if the wind is initially considerably out 
of balance the time required for the wind to adjust to 
within 5 percent of the theoretical balance velocity jis of 
the order of 10 hours under the conditions assumed. As 
an actual example, consider a wind with initial compo. 
nents uw=20 kt. and v=35 kt.; the geostrophic wind is a 
constant at v,=40 kt. and w=0. For k=0.003 mi and 
at 30° N. latitude, the wind will be in balance when 
v=34 kt. and u=15 kt. It takes 26 hours, computed by 
the methods suggested here, for the wind to adjust to 
within +5 percent of the balance values. 

Experiments were made around a hurricane using dif. 
ferent values of /. A value of 0.020 gave trajectories that 
were quite unreasonable in that they headed more directly 
across the isobars than is observed. With a value of 
0.0003 the winds acquired speeds that were much higher 
than observed and overshot the isobars leading to con- 
siderable motion toward higher pressure. A value of 
0.003 gave trajectory winds more nearly like those 
observed. 


RESULT OF FRICTICNAL ACCELERATION 


The fact that adding friction always accelerates the 
air parcel toward a certain balance value leads to an 
interesting corollary. This is that if several winds of 
different values (both in speed and direction) are all 
subjected to the same geostrophic field and this field is 
changing with time the different winds will approach the 
same value because of frictional acceleration. This is 
not necessarily the value where balance requirements are 
satisfied since the geostrophic wind itself may be changing. 

The results of a test of this corollary are shown graph- 
ically in figure 8. The four initial winds, A, B, C, and D, 
are shown at the bottom left of the figure. After 6 hours 
in a geostrophic wind field which appears as at upper 
center, their values became as at the lower center, and 
after the next 6 hours with geostrophic wind as at upper 
right the final values are as at lower right. Note that at 
the end of each of the two 6-hour intervals the differences 
among the several winds are progressively reduced. 

Trajectories in several hurricanes were computed by this 
method, using actual data for pressure fields. The results 
are encouraging. The use of /=0.003 mi.-? yielded tra 
jectories that in most respects are similar to trajectories 
drawn from actual wind data where available. These 
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GEOSTROPHIC WIND 
DURING SECOND 6 HOURS 





WINDS AFTER 6 HOURS 
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Ficure 8.—Illustration of tendency of different initial winds subjected to the same changing pressure field to approach the same 
value because of frictional acceleration. 


trajectories together with evaluations of & itself will be 
discussed at greater length in a subsequent article. 
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THE WEATHER AND CIRCULATION OF MAY 1959 
INCLUDING AN ANALYSIS OF PRECIPITATION IN RELATION TO VERTICAL MOTION 


CHARLES M. WOFFINDEN 
Extended Forecast Section, U.S. Weather Bureau, Washington, D.C. 


]. INTRODUCTION 


The May pattern at the 700-mb. level over the United 
States was one of large amplitude with a trough in the 
West and ridge in the East. Except for a short break, 
this system persisted throughout the month, bringing 
repeated outbreaks of cold maritime air into the western 
states and unusual warmth in the East. The well-marked 
inean frontal boundary thus established gave rise to fre- 
quent storminess and abundant rainfall over much of the 
country. 


2. MEAN CIRCULATION AT 700-MB. 


Perhaps the most striking aspect of the 700-mb. mean 
map for the month (fig. 1) was a tendency toward blocking 
activity. One blocking High was centered over north- 
eastern Siberia and another northeast of the Caspian Sea. 
In addition, large amplitude ridges with sizeable positive 
height anomaly centers were located over Alaska and 
Britain. Since both these features were roughly out of 
phase with the lower latitude circulation, the pattern 
which developed was strongly difluent over the western 
coast of North America and Europe and confluent down- 
stream over eastern North America and Eurasia. The 
Pacific westerlies to the south of the block were relatively 
flat with a tendency toward a weak trough off Japan and 
another in the central Pacific. This latter feature, though 
weak in the mean and variable during the month, was in 
a position to support the trough in the western United 
States and its disappearance about mid-month led to the 
only major (though transitory) change from the domi- 
nant weather pattern for the month. 

The pattern depicted in figure 1 differed in several 
significant aspects from that of April [1]. The change 
chart (fig. 2) indicates a marked increase in blocking 
activity with large rise centers over northern Asia and 
the eastern Atlantic. Also, the deepening in the western 
United States together with the ridging in the East were 
sufficient to reverse the circulation regime over the country. 

Since blocking played a major role during the period, 
its effect was studied by preparing harmonic analyses for 
the meridional component of the 700-mb. flow around 
each latitude circle. The graphs of amplitude versus wave 


number at 60° N. and 40° N. were computed following the 
procedure described in [4] and are displayed as figure 3. 
In a characteristic blocking situation with large features 
at high latitudes one would expect amplitudes in the low 
wave numbers to be dominant. On the other hand, the 
relatively short wave lengths at lower latitudes would be 
expected to favor higher wave numbers. To a degree this 
is substantiated by the data of figure 3. At 60° N., for 
example, the amplitude of the 2nd harmonic exceeds that 
calculated from the normal map [8], while all others, 
except wave 5, are less than normal. 

At 40° N. on the other hand, the amplitude of harmonies 
5, 6, and 7 were clearly greater than normal, though 1, 2, 
and 3 still retained values about average in magnitude. 
The curves at 40° N. display the double maximum de- 
scribed by several authors as characteristic (c.f., [9], [10]) 
with a primary maximum in the longer wave lengths and 
a secondary maximum near wave 6 or 7. This distribu- 
tion is also discernible in the curve representing the normal 
but with greatly diminished amplitude in the higher wave 
numbers. The tendency for the shorter wave lengths to 
predominate is further illustrated by comparing the per- 
centage of the total variation in the pattern [4] accounted 
for by wave 5, 6, and 7 in May 1959 (26 percent) with the 
corresponding percentage for the normal (4 percent). 
Thus the waves within this band were unusually energetic 
this month, and even wave 7 was sufficiently strong and 
persistent to account for a surprising 12 percent of the 
total variation in the pattern. 

The zonal index underwent a double oscillation during 
the month (fig. 4). The first was minor in nature, but 
the second was more pronounced and amounted to a major 
index cycle for this season of the year. It began when the 
zonal index briefly reached 10 m.p.s. for the 5-day period 
ending May 11, dropped to a low value of 4.4 m.p.s. for 
the period ending the 27th, and then recovered rapidly, 
again exceeding 10 m.p.s. by June 6. 


3. MEAN TEMPERATURE IN THE UNITED STATES 


The change in circulation over the United States from 
April to May brought about a corresponding reversal mn 
the temperature pattern. The observed temperature dis: 
tribution in May (fig. 5) was one of sharp contrasts with 
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Fievre 1.—Mean 700-mb. height contours (solid) and departures from normal (dotted) (both in tens of feet) for May 1959. Over 
the United States the pattern of trough in the West and ridge in the East represented a reversal from the preceding month. 
Blocking was active, as indicated by the extensive areas of positive anomaly over polar and subpolar regions. 


cold maritime air masses predominating in the mean 
trough in the West and warm tropical air in the eastern 
ridge. In the West this resulted in a regime which was 
cooler in May than in April by two classes over a broad 
irea covering most of the western Plateau (fig. 6). At 
Fresno, Calif., the temperature dropped 4 classes from 
much above normal in April to much below in May. 
Maximiimn warming, on the other hand, occurred in the 
eastern ‘lexas-Louisiana area, with a finger of two-class 
changes extending northward to the Lakes. 


At Ely, Nev., the average monthly temperature of 47.4° 
F. was the coldest of record for May. On the other hand, 
Muskegon, Mich., experienced the warmest May since 
records began with an average 60.5° for the month. 


4. EVOLUTION OF BLOCKING 


The tendency for blocking to be frequent in May has 
been discussed by several authors [5, 7], and this month 
was no exception. Two remarkable surges occurred, one 
of which influenced primarily the eastern sector of the 
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igure 2.—Change of mean 700-mb. height departure from normal 
(tens of feet) from April to May 1959. Blocking increased over 
the Siberia-Alaskan region and also over the eastern Atlantic. 
The falls in the western and rises in the eastern United States 
resulted in a change of phase from the April pattern. 


Northern Hemisphere and the other the western. To fol- 
low their evolution, the series of 5-day mean maps ap- 
proximately a week apart (fig. 7) has been prepared. The 
history of the first blocking wave began with the subpolar 
High over Alaska on the first chart (fig. 7a) of this series. 
This center shifted rapidly to northeastern Siberia on the 
subsequent map (fig. 7b), and its influence can be traced 
even farther westward as a band of positive height anom- 
aly extended across the Arctic basin and joined the next 
upstream ridge near Novaya Zemlya. Thereafter, the 
original center sank southward to the neighborhood of 
Kamchatka (fig. 7c) and then appeared to join with the 
subtropical ridge in the western Pacific (fig. 7d). 
However, the blocking surge, as distinguished from this 
center, can be traced on to the west as heights subsequently 
rose over 600 ft. in the Scandinavian region from the 
second to the third chart of the series, resulting in a 
marked retrogression of the Russian blocking ridge to that 


region. The North Atlantic was next affected, and a 


closed anticyclone appeared just north of Britain (fig. 7d). 
By the end of the month, this surge seemed to come to an 
abrupt end as heights fell rapidly at northerly latitudes 
from Europe to North America. 

So far this phenomenon has been treated in terms of 
ridges, but one could equally well direct attention to the 
remarkable sequence of Lows cut off or depressed to lower 























WAVE NUMBER 


Figure 3.—Amplitude in tens of feet of the zonal harmonic wave 
components at 60° N. and 40° N. for the 700-mb. mean chart for 
May 1959. At 60° the maximum amplitude was in the second 
harmonic. At 40° N. the graph exhibits a double peak with 
maxima in waves 2 and 7. The fifth, sixth, and seventh har- 
monics were markedly above normal (dashed), illustrating the 
tendency for short wave lengths to be associated with blocking 


latitudes: first over the sea of Okhotsk (fig. 7a), then just 
north of Lake Baikal and over the Black Sea (fig. 7b), 
then near the Adriatic and Central Atlantic (fig. 7c), and 
finally over the Bay of Biscay and the Labrador Shelf 
(fig. 7d). This long list illustrates the widespread in- 
fluence which such a blocking impulse can exert and high- 
lights its possible importance as a tool in longer period 
forecasting. From its position over Alaska at the be- 
ginning to that north of Scotland near the end of the 
period it traversed approximately 205° of longitude in 1! 
days or about 11° of longitude per day, an unusually rapid 
pace. 

For the western sector of the hemisphere, however, the 
more important blocking complex was that centered off the 
Labrador coast as the month opened (fig. 7a). In succes 
sive weeks this surge, spreading upstream : 

1. Amplified the United States ridge and cut off the 
Low near Newfoundland (fig. 7b). 

2. Amplified the ridge in western Canada and cuit off the 
Low near James Bay (fig. 7c). 
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FieureE 4.—Time variation during May and early June 1959 of 
the 5-day mean values of the zonal westerlies in meters per 
second (plotted on the last day of the period). The zonal index 
is computed for the western sector of the hemisphere for the 
latitude belt 35°—55° N. A pronounced index cycle followed an 
earlier minor variation, with the zonal index reaching a minimum 
during the last week of May. 


3. Amplified the eastern Pacific ridge, cut off the de- 
pression in the Southwest, and began to depress the cyclone 
inthe eastern Aleutians (fig. 7d). 

4. Amplified the pattern over northeastern Siberia and 
cut off the deep system in the eastern Aleutians (fig. Te). 

This singular evolution paralleled the first such system 
inmany respects and can be described as typical of many 
blocking sequences. It was intimately connected with the 
index cycle previously described, the low point of the latter 
being reached as the blocking reached maximum coverage 
ind development in the western sector of the hemisphere. 
Also, subsequent continued retrogression of blocking into 
Asia in early June allowed the zonal index to again climb 
to above 10 m.p.s. and complete the cycle. The average 
speed of retrogression of blocking in this case was about 
\* of longitude per day, which is considerably slower than 
the first case described. 


5. WEEKLY WEATHER AND CIRCULATION OVER 
THE UNITED STATES 


The wave train over the United States was progressive 
for the first portion of the month. The trough initially 














Figure 5.—Observed temperature anomaly for May 1959 expressed 
in five classes defined so that near, below, and above normal 
each comprises one-fourth of the monthly mean temperatures for 
May over the 30-year period 1921 through 1950. The remaining 
fourth is divided equally between the much above and much 
below categories. A strong east-west temperature gradient 
existed, with very cold conditions prevailing over the West 
(with the exception of coastal California) and very warm 
weatber over roughly the eastern two-thirds of the country. A 
marked frontal system separated these extremes, and wave 
activity resulted in widespread precipitation and severe storms. 
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Figure 6.—The number of classes the anomaly of temperature 
changed from April to May 1959. Marked cooling took place 
over a broad area of the West, with warming over most of the 
remainder of the United States except for the northeastern 
States. 


along the west coast (fig. 7a) marched steadily eastward 
and by the fourth map of the series (fig. 7d) was well out 
in the Atlantic. It has been mentioned that blocking 
served to amplify first the ridge in advance of this trough 
(fig. 7b) (thus advecting warm moist air into the Missis- 
sippi Valley) and subsequently the Canadian ridge in its 
rear (fig. 7c) (thus bringing down unusually cold air). 
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FicureE 7.—Five-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for (a) 
April 30-May 4, (b) May 5-9, (c) May 12-16, (d) May 19-23, 

3B. 0+ and (e) May 28-June 1, 1959. A mean trough, which entered 
the west coast early in the month, marched regularly across the 
country and passed through the East about mid-month. There 
after it progressed rapidly into the Atlantic, and the pattern 
over the United States reverted to the original one with a deep 
trough in the West and ridge in the East. The May index 
cycle appears related to the blocking wave which emanated from 
the Atlantic toward the end of April and can be traced westward 
to the Pacific by late May. 
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FieurE 8.—Departure of average temperature (° F.) from normal 
for the weeks ending at midnight local time (a) May 3, (b) May 
10, (c) May 17, (d) May 24, and (e) May 31, 1959. The cold 
air which entered the Pacific States early in May marched 
steadily eastward until it overlay most of the East by the week 
ending May 17. By that time, however, above normal tempera- 
tures appeared in the Southwest and subsequently overspread 
the East by the following week (ending May 10). Thereafter, 
this regime became very persistent and largely determined the 
character of the temperature anomaly distribution for the 
month. (From Weekly Weather and Crop Bulletin, National 
Summary, vol. XLVI, Nos. 18 through 22, 1959.) 


Along the boundary between these contrasting air masses 
severe frontal weather developed with thunderstorms, 
hail, high winds, and a number of tornadoes. As a result 
precipitation was heavy over much of the area between 
the Rockies and the Alleghenies, and was locally excessive 
over the Mississippi and Ohio Valleys. The eastward 
progress of this frontal boundary can be traced with the 
aid of figure 8 which shows the cold air reaching the East 
luring the week ending May 17 (fig. 8c). This brought 
the only cold snap of the month to that region. There- 
after tle basic circulation pattern of the month (fig. 1) 
became reestablished as a deep trough entered the West 
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Bannd om pectmnary ‘einwephic rope 
(fig. 7d) and the ridge over the Rockies sheared from its 
previous connection with the blocking ridge over north- 
western Canada (fig. 7c) and amplified as it moved 
rapidly to the east coast (fig. 7d). It is noteworthy that 
the East experienced a complete reversal from very warm 
the second week (fig. 8b), to very cold the next (fig. 8c), 
and back to the original regime the fourth week (fig. 8d). 
Numerous records of daily maximum and minimum 
temperatures were set, including maxima during the early 
period of 99° F. at Huron, S. Dak., on the Ist; 95° 
Louisville, Ky., on the 4th; 95° at Augusta, Ga., on the 
3d; and 82° at Portland, Maine, on the 11th. These were 
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followed by record daily minima of 24° at Huron on the 
15th; 38° at Wilmington, Del., on the 16th; 46° at 
Augusta, Ga., on the 16th; and 31° at Portland, Maine, on 
the 18th. The warm period which followed was not as 
record-breaking, though Portland, Maine, did manage to 
reach a new record maximum for the 28th of 88° F. 


6. WEEKLY HARMONIC ANALYSIS 


A harmonic analysis was also prepared for each of the 
charts in figure 7, and the results for 60° N. and 40° N. are 
presented in figure 9. As previously discussed, it was ex- 
pected that blocking would tend to accentuate the longer 
wave lengths at high latitudes and shorter wave lengths at 
lower latitudes. However, the hemispheric circulation 
pattern is often complicated by a tendency for low zonal 
index in one area to be offset by high index elsewhere [2]. 
This suggests that. perhaps fairly even distribution over a 
wide wave length spectrum should be more characteristic, 
particularly at lower latitudes. In a general way this is 
substantiated by the graphs at 40° N. For example, dur- 
ing the low index period May 19-23, the amplitudes are 
quite regularly distributed with an amplitude of 50 ft. in 
the 9th harmonic, a sizeable figure for that wave number. 
Apart from the relatively evenly distributed amplitudes at 
40° N. over a rather wide frequency range, the maxima 
tended to cluster about waves 1, 2 or 3, 6 or 7, and (in the 
case of lowest zonal index) 8 or 9. In fact, these latter 

rave numbers, which correspond to wave lengths of only 
45° and 40°, respectively, are quite characteristic of those 
regions where blocking was most active as, for example, 
over the western Atlantic, May 12-16 (fig. 7c) ; the eastern 
Pacific, May 19-23 (fig. 7d), and the western Pacific, May 
28 to June 1 (fig. Te). It is interesting that wave lengths 
of this magnitude, which are of the dimensions of daily 
systems, at times prevail even in the mean. It is likely 
that they occur mainly during periods of low index. 

The effect of blocking at high latitudes, however, does 
seem to be to accentuate amplitudes in the low wave num- 
bers. Wave 1 appears to be favored if only one blocking 
system exists or is strongly dominant. This was the case 
during the period April 30-May 4 when the Eurasian 
High was much stronger along the 60th parallel than its 
counterparts over Alaska and south of Greenland. How- 
ever, the latter two did make substantial contributions to 
the amplitude of waves 2 and 3 which were also unusually 
large during this period. 

On the other hand, when dual blocking Highs exist 
across the pole from each other, and this is often the case 
during periods of lowest index, wave 2 is principally 
affected. The graph at the 60° parallel for May 19-23 
(fig. 9) may be cited as an example. 


7. PRECIPITATION 


A large portion of the United States enjoyed adequate 
rainfall during May (fig. 10). Heaviest amounts fell over 
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Figure 9.—Amplitude in tens of feet of the zonal harmonic wave 
components at 60° N. and 40° N. for 5-day mean 700-mb. charts 
for the periods indicated. At 60° N. the low wave numbers 
dominated much more strongly than at 40° N. On the other 
hand, amplitudes at 40° N. were much more evenly distributed 
over a fairly wide spectrum, with some tendency for maxima 
to appear at the very long wave end and again in the wave 
number band 5 through 7. 





Percentage of Norma! Precipitation for May 1959 





Based on preliminary telegraphic reports 








Figure 10.—Percentage of normal precipitation for May 1959. 
Precipitation was heavy over much of the Nation. Note the 
continued dryness in the Southwest and near the border of 
Montana and Wyoming with the Dakotas. (From Weekly 
Weather and Crop Bulletin, National Summary, vol. XLVI, No 
23, 1959.) 


the Central Plains in advance of the mean trough position 
for the month (fig. 1), in the Northwest, and also in the 
Southeast. The heavy precipitation in the Sout!iast 
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FicurE 11.—Estimate of 30-day mean vertical motion for May 1959 in units of mm./sec., with positive values representing upward 
motion, obtained by: (a) advecting the mean 850-500-mb. thickness with the mean 500-mb. flow, (b) averaging 60 twice-daily values 
at 600 mb. computed from initial data by the JNWP baroclinic model, (c) applying the JNWP baroclinic model to the observed 
mean monthly 500-mb. contours and 850-500-mb. thickness data, (d) applying an equivalent barotropic model to the mean 
500-mb. contours for the month. Only positive values of vertical motion were included. 


curred during the latter half month, partly in connection 
with frequent locally heavy showers and thunderstorms, 
and also with a large contribution from tropical storm 
Arlene, particularly along the coastal section of Louisiana. 
Total rainfall at Rome, Ga., for example, reached 9.54 
inches to set a new record since 1856. 

The continued dry weather in the Southwest aggravated 
the existing drought situation in that area. Tucson, Ariz., 
experienced its first rainless May since 1949. The 0.32 in. 
recorded from January 1 to May 31 made this the driest 
comparable period since records began in 1868. Also, 
Oakland and Sacramento, Calif., with only traces, re- 
ported this month to be among the driest of record. Since 
January 1 El] Paso, Texas, experienced only 0.73 in. (43 
percent of normal) while Wichita Falls, Texas, reported 
this May to be the 10th consecutive month with subnormal 
precipitation. 

New F ngland and the Central Atlantic States were also 
relatively dry under the influence of anticyclonic flow and 


offshore winds. Hartford, Conn., and Scranton, Pa., for 
example, reported the driest May of record with monthly 
totals of only 0.73 and 0.77 in. respectively. 


8. SPECIFICATION OF PRECIPITATION BY VERTICAL 
MOTION 


The advent of the electronic computing machine has 
made possible routine computation of vertical motion in 
the middle troposphere. Several authors have recently 
discussed the relationship between this parameter and ob- 
served rainfall. Sanders [6], for example, investigated 
the occurrence of cloudiness and precipitation at certain 
cities in the United States and discussed their relation- 
ship to concurrent values of the dynamical, orographic, 
and frictional contributions to the vertical motion. 

In the following experiment, several estimates of ver- 
tical motion were computed (fig. 11) and compared to 
the observed rainfall anomaly pattern (fig. 10). Since 
production of precipitation is a discontinuous and com- 
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plex process dependent upon a number of parameters, use 
of this single variable can hardly be expected to account 
for the whole pattern. 

The first of these charts (fig. 1la) depicts the advection 
of the mean thickness field (850 to 500 mb.) by the mean 
500-mb. flow for the month. By assuming the thickness 
field to be stationary, which is a reasonable assumption on 
the monthly mean basis, the result can be interpreted in 
terms of upward motion. Comparison with the observed 
rainfall distribution (fig. 10) reveals a reasonable simi- 
larity between the two patterns. The relatively wet area 
in the central portion of the country appears to be fairly 
well delineated from the drier northeastern states by the 
line of neutral advection. In the Southeast, however, the 
advective pattern, though positive, does not support as 
much rainfall as occurred, particularly in eastern Georgia. 
In the mountainous areas of the West where orographic 
effects are likely to dominate, one is not surprised to find 
a poor relationship. In this region, however, although 
precipitation was greater than normal for the month, it 
totalled only one to two inches in absolute amount.* 

The second chart (fig. 11b) was obtained by averaging 
the twice daily values of observed vertical motion at 600 
mb, as computed from the two-level baroclinic model cur- 
rently used by the Joint Numerical Weather Prediction 
Unit (JNWP). Also in this case, the resultant pattern 
bore a close correspondence to the observed rainfall pat- 
tern. The band of heavy precipitation extending from 
western Texas to the Great Lakes was associated with 
average upward motion as were the heavy rains in the 
Southeast. On the other hand, average descent in the 
Northeast and extreme Southwest was apparently quite 
effective in keeping those regions dry. As in the previous 
case, the correspondence was poor over the Northwest. 

Figure 11c shows the vertical velocity obtained by using 
the monthly mean height at 500 mb. and thickness from 
850 to 500 mb. as input data to the JNWP baroclinic 
model and running for one iteration. Thus this chart is 
closely akin to the previous one and the two would be the 
same if relationships between circulation and vertical 
motion were linear. The two charts are indeed quite 
similar, and most of the comments of the preceding para- 
graph apply in this instance also. The area of largest 
discrepancy over the United States lay over the States 
bordering Lake Michigan, where values computed in this 
manner were larger by 2 or 3 mm./sec. It is of interest 
that in both cases mean troughs were associated with rising 
motion in advance and sinking motion behind the trough 
line. 

Figure 11d has been included to test the extent to which 
a simple model can capture the overall aspects of the 
monthly precipitation field, given the flow pattern. In 
this case the barotropic model described by Namias and 


*For total precipitation in inches see Weekly Weather and Crop Bul- 
letin, vol. XLVI, No. 22, June 1, 1959. 
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Figure 12.—Number of days in May 1959 with fronts of any type 
within unit squares (with sides approximately 500 miles). All 
frontal positions are taken from Daily Weather Map, 1300 psy. 
Areas with 15 days or more with fronts are stippled. Active 
fronts were frequent over the central portion of the country. 


collaborators [3] was applied to the mean 500-mb. flow 
except that the basic current was omitted (see appendix 
for method of computation of vertical motion). The 
model was designed to include positive vertical velocity 
only, on the assumption that precipitation might be better 
related to total upward motion than to the arithmetic 
average. Figure 1ld was the result after one time step. 
Though the pattern was ill defined because the values were 
small, it bore a distinct resemblance to figure 11b, with 
cells of maximum positive vertical motion in the South- 
west, the central Gulf States, and western Great Lakes 
Region. Each of these cells covered a region of consider- 
able precipitation, though the area of very heavy rainfall 
along the Nebraska-Iowa border was not well indicated. 
Since the chart series, figure 11, constitutes but a single 
case, little in the way of conclusions can be drawn. How- 
ever, they do suggest that continued research aimed at 
relating mean precipitation to a vertical motion parameter 
might be fruitful. 

In general the charts (fig. 11) define an axis of maxt- 
mum ascending motion extending roughly from western 
Texas to the Great Lakes. It is interesting that this 
agrees remarkably well with the position of the most fre- 
quent frontal activity during the month (fig. 12). 


APPENDIX 


The vertical motion w was estimated from the following 
equations derived by P. F. Clapp (personal communica 
tion). 

The equivalent barotropic atmosphere is characterized 
by invariant wind direction with height so that if V: 





-—s = wm Co - —_— 








1969 


SA 


v 


~~ OF, 


é 


> 





low 
dix 
The 
city 
tter 
etic 
tep. 
ere 
rith 
ith- 
kes 
ler- 
fall 


gle 
ow- 
at 
ter 


eri 
his 


ing 
ca- 


IS 


May 196° 


the wir 1 at one level, then the wind at a different level is 
v=AV., where A is a function of height. Thus, if the 
subscript 1 is taken to refer to 500 mb. and the subscript 
9to 75 mb. and we assume conditions at 500 mb. to be 
nondivergent, the simplified vorticity equation at 500 and 
750 mb. can then be written: 


OF EV Vt Mi VI=0 (1) 
Of, o\_ 
Bt Va: Vist We- WIS ($2) =o (2) 


where the symbols have their usual meaning and a is 
the vertical velocity in terms of pressure coordinates. 
Combining (1) and (2) and noting that ¢,=A¢,, we obtain 


w 1 ; K : 
(So) =p ADV =5 V,-Ve, (3) 


where K is a constant. 


If we approximate ($<) by ong and assume that w at 
2 1 


1000 mb. (denoted by ws) is zero, then 
w= —9 KV,-Ve, (4) 


From the hydrostatic equation 


dp 
dt 


w= 


) F=—pPigwi (5) 
1 


Thus w, can be estimated from (4) and (5) subject to the 
shortcomings of the equivalent barotropic model and the 
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assumptions made above. In figure 11d it has been scaled 
in units of millimeters per second. 


3. 


10. 
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Publication of Weather Notes 


Many years ago the Monthly Weather Review published detailed eyewit- 
ness accounts of exceptional storms. These accounts both enrich the meteor- 
ologist’s knowledge of storms and provide him with particular details that 
cannot be found elsewhere. Because such information bears directly upon 
questions the meteorologist must attempt to answer about weather phenomena 
(for example, the identification of storms as tornadoes), and because the 
information has potential value in both the research and service programs of 
the Weather Bureau, publication of eyewitness accounts and brief analyses 
of exceptional storms and other meteorological phenomena was resumed in 
the April 1955 issue. They appear from time to time under the heading 
“Weather Notes.” 

Contributions to these “Notes” are invited from readers of the Review. 
There is no limitation placed on length of description but it is expected that 
most will be short accounts. Any weather peculiarities, whether storms or 
other phenomena, are acceptable subject matter. Material should be addressed 
to Editor, Monthly Weather Review, U.S. Weather Bureau, Washington 25, 
D.C. 
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